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Bei v/w = 1/2 werden diese Kineton Energien der Abgase beitler Fahrzeuge 
gerade gleich, wenn auch ihre Géschwindigkeitsrichtungen entgegengesetzt sind. 
Bei, grosseren Fahrgeschwindigkeiten ist die Rakete aus denselben Griinden giin- 
stiger und wird es fiir immer héhere v/w inimer mehr, da die Wirkungsgrade 
des Raketen- und des Radantriebes sich wie 2 v/w verhalten. . 

Dies deshalb, weil mit v/w > oo bei unter sich gleich werdender kinetischer 
Energie der Abgase beider Antriebe® die Schubleistung des Raketenantriebes 


2 3 
iw —v/V2gH 
Figur 1 Figur 2 
Grundgleichungen des Radantriebes Aussere Wirkungsgrade des Radantriebes 
und des Raketenantriebes. und des Raketenantriebes. 


wv dm,/dt immer weiter wachst, jene w?/2- dm,/dt des Radantriebes aber kon- 
stant bleibt. 

Selbst wenn Schienen durch den Weltraum fiihrten, ware es demnach meist 
wirtschaftlicher, dort nicht mit radangetriebenen, sondern mit raketenangetrie- 
benen Fahrzeugen zu verkehren. 


Summary 


It is a general believing that propulsive efficiency of wheel-propulsion for 
vehicles is always unity, and thus superior to propulsive efficiency of rocket-pro- 
pulsion. By some elementary mechanical considerations can be shown that this 
is only true for slow drive velocities. 

With growing velocities, propulsive thrust of wheel-propulsion decreases, 
while it remains constant for rocket-propulsion; propulsive performance remains 
constant for wheel-propulsion, whilst it grows for rocket-propulsion with increas- 


ing velocity. 


As a consequence thereof, propulsive efficiency of wheel-propulsion decreases 
much faster than that of rocket-propulsion with growing velocity, so that efficiency 
of rocket-propulsion is far higher at high velocities. 


(Eingegangen: 8. Oktober 1953.) 
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An Arithmetic Unit for Automatic Digital Computers 


By Joun R. Stock, Ziirich’) 


A number of moderate-speed, automatic digital computers using a magnetic 
drum for storage have been built or designed?). In this communication, certain 
general requirements on the arithmetic unit of such a computer are presented, as 
have been stipulated for the new computer at the Swiss Federal Institute of 
Technology?). Then, to fulfill these requirements, a suitable set of arithmetic 
operations is proposed and a design of the arithmetic unit is sketched. 

Some scheme should be employed in the computer to permit the represen- 
tation of a broader range of numbers than is feasible with a fixed decimal point. 
When a wide variety of practical problems is likely to be encountered, this is a — 
matter of great convenience for the mathematician using the machine. The float- 
ing decimal point is very suitable; it has the particular advantage of being — 
thoroughly familiar. Other schemes, such as the use of fixed decimal point num- 
bers | *|< 1 together with their reciprocals, either complicate the arithmetic 
operations excessively or require similar shifting operations to those which char- 
acterize computation with a floating decimal point. 

Computation should be done in the decimal system in preference to the binary. 
With a floating point, automatic conversion to and from the binary system is 
somewhat difficult, and the equipment used for this purpose might just as well 
be used for computation directly in a coded decimal system. The decimal system 
also facilitates the mathematical investigation of errors. 

It must be possible to calculate with double or multiple precision. This feature 
is perhaps more a matter of principle than of absolute necessity, but lack of it 
would be deemed a definite deficiency. Certain problems arise, however, in com- 
bining double precision with a floating decimal point, since the former requires 
the retention of information which tends to be dropped by the latter. Besides 
operation with double precision, the arithmetic unit must be able to carry out 
calculations upon commands. In particular, this implies being able to separate 
a given command into its component parts (operation and address), to transform 
these by addition or substitution, and to reassemble a new command. 

In any automatic computer, the overall speed of calculation is limited by the 
average access time to storage, which for a magnetic drum is of the order of 
5 to 10 milliseconds. There is little to be won by building an arithmetic unit which 
is much faster than this, unless the mathematician is willing to do minimum access 
programing. On the other hand, little is lost if less frequent operations, such as 
multiplication and division, require several times the access time. This fact speaks 
in favor of a floating decimal point, since the floating decimal point tends parti- 
cularly to slow down the operation of addition, which normally is disproportio- 
nately fast. Thus the machine accomplishes more without serious loss of overall 
speed. A programed floating decimal point, however, is much slower and would 
hardly be permissible in a moderate-speed machine. Again, computation directly 
in the decimal system is preferable to programed conversion to the binary system. 


1) Institut fiir angewandte Mathematik der ETH. 
2) Papers describing several of these are listed in the Bibliography, page 171. 


3) A. P. SPEIsSER, Projekt einer elektronischen Rechenmaschine am Institut fi angewandte 
Mathematik der ETH., ZAMP 4, 317 (1958). 
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_ By a similar line of thought, all routine modification of the address of a command 
should be accomplished in the sequence unit under the control of some of the 
digits of the operation part of the command. Computation in the arithmetic unit 
upon commands should be restricted to special cases?). 

Last but not least, the operations which the machine can perform should be 
in such a form that they can be readily understood and learned by any engineer 
or mathematician who may have an opportunity to use the machine. 

To realize these requirements, it is proposed that three types of ‘words’ be 
utilized, each consisting of 15 decimal digits plus one digit reserved for checking 
purposes. The three types are 
(1) A number with floating decimal point, ax 108. Range of exponent b: + 99. 

Range of mantissa a: + 0-99999999999, 

(2) Anumber with fixed decimal point. Range of mantissa: + 0-99999999999999, 
(3) Two one-address commands, each with three-digit operation and four-digit 
address. 

In general, an arithmetic operation combines two words, one of which is 
taken from a register MR in the arithmetic unit and the other from storage. The 
result of the operation is usually only one word, which remains in two registers, 
MR and AC. Occasionally the result may be two words, in which case one remains 
in MR and the other in AC. In the following list of the more essential arithmetic 
operations, the letter 7 indicates that the address of the command is utilized, and 
(7) represents the number in storage position 7. 

Operations defined for both numbers with fixed and floating decimal points: 

A n: Clear MR and AC, and insert (7) into them. 

Sn: Store the number in MR into position », and then interchange the 
numbers in MR and AC. This interchange makes the number in AC 
available for the next operation should MR and AC contain different 
numbers. 

—1: Multiply the number in MR by —1; insert the result as usual into MR 
and AC. 

|«|: Take absolute value of number in MR. 

Operations defined only for numbers with a floating decimal point: 
N: Normalize the number in MR; that is, shift mantissa left with corre- 
sponding change of exponent until most significant digit is not zero. 

+ n: Add (x) to number in MR with round-off. 

— mn: Subtract with round-off. 

x nm: Multiply with round-off. 

+ mn: Divide with round-off. 

Operations defined only for numbers with a fixed decimal point: 
+ -- »: Add; stop if sum greater than 1. 
— — nv: Subtract. 
x x n: Multiply. Of the double length product, the more significant half 
remains in MR;; the less significant, in AC. 
+ + mn: Divide. Stop if divisor is less than dividend. Quotient remains in MR; 
remainder in AC. 

L n: Left shift of the 14 significant digits by places; that is, multiply by 
10” and drop any overflow = 1. 

Rm: Corresponding right shift. 


1) H. RuTIsHAusER, Automatische Rechenplanfertigung bet programmegesteuerten Rechenma- 
schinen, Mitteilungen aus dem Institut fiir angewandte Mathematik der ETH., Nr. 3 (Birkhauser, 


Basel, 1952). 
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The various special cases, such as overflow of an exponent, normalization of 
zero, and division by zero, will not be discussed here. The arithmetic unit must 
also supply certain information to the sequence unit for conditional calls. Useful 
are the algebraic signs, mantissa identical to zero, and threatening overflow, 
that is, exponent = 50 or mantissa => 0-5 for floating and fixed decimal points 
respectively. 

The inclusion of operations on fixed decimal point numbers makes possible 
computation with multiple precision and computation with commands. In parti- 
cular, the left and right shift operations allow the extraction of any digit or 
group of digits in a word. There is the interesting possibility of using the left 
from storage 
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Figure 1 
Block Diagram of Arithmetic Unit 


shift to clear the exponent out of a number, leaving only the mantissa as a 14- 
digit number with a fixed decimal point. This is useful in the computation of the 
elementary functions. 

The proposed arithmetic unit to carry out these operations is basically of 
conventional design. It contains three circulating registers, known as the multi- 
plicand, the accumulator, and the multiplier, each holding one word of normal 
length. Numbers coming from storage enter the MD register, and those being 
stored are taken from MR. The coded decimal digits are treated in parallel, but 
the 16 digits in a word are in serial order. The more important logical circuits 
are the main adder (adds two decimal digits plus carry), the decimal-binary 
converter (recodes a two-digit decimal number into a hexadecimal digit, setting 
numbers larger than 15 equal to 15), and the small binary adder (subtracts a 
unit from a hexadecimal digit.) In the block diagram shown, the contents of the 
three registers flow from left to right, closing the loop over the edge of the 
figure. Circles represent logical circuits; squares, unit delays. 

The arithmetic unit differs from conventional designs, though, in several 
ways. The exponents are not handled in separate registers from the mantissas, 
but in the same ones at different times. This represents a saving of material and 
permits the use of the exponent places for mantissa digits in numbers with a 
fixed decimal point. It has been found not necessary to include empty digits in 
the numbers for computational purposes, since sufficient space is available in 
the decimal digit used to represent the algebraic sign. To control the number 


SE 
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of repeated additions required by a multiplier digit during multiplication, no 
counter is employed, but instead the MR digit itself is reduced by one for each 
addition. A similar method is used in shifting operations, in that the number of 
shifts is inserted as a hexadecimal digit into the sign digit of the MR register. 
For each shift of one place, the hexadecimal digit is reduced by one. This permits 
shifts up to fifteen places, which is actually more than necessary. Finally, there is 
no counter to indicate the end of a multiplication according to the number of 
digits in the multiplier. For this purpose a code combination which does not 
represent a decimal is inserted into the MR sign digit and is simply shifted with 
the other digits until it reappears. In this way the same control may be used for 
the 11-digit multiplication x and the 14-digit multiplication x x. 

In more detail, the operation + proceeds in the following manner. First AC 
and MD are cleared and the addend is inserted from storage into MD. The expo- 
nents are subtracted in the AC register, the absolute value of the difference 
being recoded as a hexadecimal and inserted into the MR sign. Depending on the 
sign of the exponent difference, the proper mantissa is brought into AC and 
shifted to the right. The addition proper may then take place, followed by round- 
ing off, taking care of any overflow, and inserting the result into MR and AC. 
Multiplication and division are similar. For the multiplication proper, the simple 
method of repeated additions is somewhat too slow. A substantial increase in 
speed is gained by subtracting instead of adding for MR digits = 5 and then 
adding one extra time for the next digit!). For example, the multiplier 137 is 
treated as 1 4 —3. For this purpose the multiplier is recoded by allowing it to 
circulate once through a special circuit. Division is done by repeated subtraction 
using a small circuit to indicate when the remainder becomes smaller than the 
divisor. The quotient is built up in the MR register as straight binary coded 
decimal digits, which may be readily converted to the 2* code used otherwise 
throughout the machine by being allowed to circulate once through the main 
adder. 

Automatic checking is included. Associated with each number is a check digit 
congruent to the number modulo three. This is the simplest check which has the 
properties of being dependent upon all the digits of the number and of being 
arithmetically invariant. The latter property allows the check digit to circulate 
through the main adder and permits arithmetic operations to be carried out on it. 

With the relatively low frequency of 32 kc, the circulation time of a word 
is 0-5 ms. This gives average computing times of 5 ms for +, 20 ms for x and x x, 
30 ms for + and ~~, and 1 or 2 ms for the other operations. These times are 
in good accord with the access time of a magnetic drum. 
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Zusammenfassung 


Fiir eine programmgesteuerte Rechenmaschine mit einer magnetischen Trom- 
mel als Speicher wird ein fiir den allgemeinen Gebrauch geeignetes Verzeichnis 
der arithmetischen Befehle gegeben und der logische Aufbau des Rechenwerks 
beschrieben. Vorgesehen ist das Rechnen sowohl mit gleitendem Komma (Zahlen 
der Form a- 10°) als mit festem Komma, was das automatische Berechnen von 
Rechenplanen sowie das Rechnen mit verdoppelter Stellenzahl gestattet. 


(Received: October 24, 1953.) 


Uber eine Funktionalgleichung aus der Warmeleitung 


Von HERMANN WunNDT, Freiburg i. Br.?) 


Zur Bestimmung der mittleren logarithmischen Temperaturdifferenz fiir Gleich- 
und Gegenstrom-Warmeaustauscher 


At, — At, 

Mie Tt AT ATA 
wurde 1951 von W. MEHNER ein Nomogramm entworfen (BW K-Arbeitsblatt 12), 
das aus drei gleichabstaéndigen parallelen Leitern besteht, die je dieselbe Funk- 
tionsskala tragen. Werden auf den ausseren Leitern die At,- bzw. At,-Werte ein- 
gezeichnet, so findet man At,, als Schnitt der mittleren Leiter mit der Verbin- 

dungsgeraden von Af, nach Af,. 

Es taucht nun die Frage auf, welche Funktion f(x) das Verlangte leistet und 


demnach auf jeder der drei Leitern aufzutragen ist. Eine einfache geometrische 
Uberlegung zeigt, dass f(#) der Relation 


flay) = r( % — Xs _ ta) 


An (4;/%2) 


ia) (0 = 4 = ) (1) 
n (44/%9) 
zu gentigen hat?). Die reellen, stetigen Lésungen dieser Funktionalgleichung sind 
zu bestimmen. 
Vor dem eigentlichen Lésungsgang kann man sich schon die folgenden Eigen- 
schaften der gesuchten Funktion f(#) klarmachen: 
a) f(#) ist nur bis auf einen willkiirlichen Faktor A festgelegt: 


| 
2 


A f(%m) mad A 1(4) - A } (42) ; 


1) Mathematisches Institut der Universitat, Abteilung Angewandte Mathematik. 
2) Die Anregung zur Beschaftigung mit dieser Fragestellung verdanke ich den Herren Dr. K. 
Srork und Dr. W. von BULTzINGSLOWEN. 
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_ An einer nichtsingularen Stelle kann daher der Funktionswert noch beliebig vor- 


geschrieben werden. 
b) /(*) ist nur bis auf eine willkiirliche additive Konstante B bestimmt: 


[f(4) + Bl + [f(*.) + B] : 


2 


f+m) + B= 


c) Die Gleichung (1) sagt aus, dass der Funktionswert an einer Stelle x stets 
gleich dem arithmetischen Mittel der Funktionswerte an den beiden Stellen x, 
und #, ist. Es sei nun etwa ¥,< %. Man findet, dass stets 

Wy Sry S *, 
gilt, und fiir endliche *,, x, + 6 kann man sogar die Gleichheitszeichen weglassen 
(¥y ist ein echter « Mittelwert » von x, und ¥,). Daraus folgt, dass {(%) monoton ist. 
Denn ware dies nicht so, so existierten zwei Stellen v,+ ¥%, fiir die f(%,) = f(%,) 
gilt. Es gibt dann eine Stelle *y im Inneren von <x,, ¥,>, die gleichfalls den 
Funktionswert 


hat. Genau so schliesst man fiir das Intervall <%,, ¥),> usw. f(¥) kann nur dann 
konstant sein, falls man es als stetig voraussetzt. Die trivialen Falle 


{(#) = B = const 


gehéren somit zur allgemeinen Losung, die jetzt in der Form /(x) =A (x) +B 
erwartet wird. 

Fiir die nichttrivialen Falle g(x) # 0, falls es solche gibt, setzt man nun 
io Oe Onind serhalteaus (1) 992 7(0) =f) -—:7(4,),. das) heisst fir 
f(~) $ #(0) = const den Wert f(0) = +-oo. Im nichttrivialen Fall ist also f(x) an 
der Stelle 0 singular. Da aber /(*) nicht tiberall singular sein soll, so folgt, dass die 
Annahme /(%), = f(*.) fiir ¥, + ¥, zum Widerspruch fihrt. /(¥) ist monoton in 
allen Intervallen, die die Null ausschliessen. 

d) Auch die Ableitung /’(*) ist fiir nichtkonstante Losungen eigentlich mono- 
ton. Bei Beschrankung auf reelle Argumente miissen ¥, und *, gleiches Vorzeichen 
haben, da sonst der Logarithmus im Nenner von ¥y komplex wird. Es sei nun 
4, > *%, > 0; fiir ¥, und x, negativ gilt analog das Umgekehrte. Man schatzt ab 


¥4— %# *14 
eae 1 See : 
In) 45 2 
fiir alle endlichen Paare (%,, ¥,), und 
. xy aa Xo 
lim *#y= —,; 
X)¥,—> CO act 


bei endlich bleibender Differenz x, — 72. Daher sagt (1) aus, dass das arithme- 
tische Mittel zweier Funktionswerte innerhalb des zwischen diesen Stellen lie- 
genden Intervalls, von links kommend, schon vor dey Intervallmitte erreicht wird, 
Daraus folgt, dass fiir fallendes f(x) und « > 0 die die Funktion darstellende 
Kurve nach unten, fiir steigendes /(¥) dagegen nach oben gekriimmt ist. Fiir 
% < 0 gilt entsprechend das Umgekehrte. 
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Nach diesen vorbereitenden Bemerkungen kann die Funktionalgleichung nach 
folgendem allgemeinem Verfahren geschlossen gelést werden. Man setzt 


se ee =e | oe 


Flay, %q, ¥m(¥%1 %2)] = —fo(¥m) + 


wobei x, und #, als unabhangige Variable aufgefasst werden. Nun bildet man 
OF /dx, und 0F/0x, (Existenz der bendtigten Ableitungen vorausgesetzt) : 
OF Ofo 0% OF Of; OF Ofs 


he fouls Ce 
ify Oxy 0%, | Of, 0%, Of OF, 


(3) 
OFS Sdjs| (aut OF NO (OR Nia 0s | 
Of, Oxy. O40, Oe. Vein wees 
Dies ergibt im vorliegenden Falle 
fo (. et ae Bo DEN ey 
On \ In ea) ¥e Papi (eee) 2 ak, ; (4) 
plies | : sree eee Wey 
| Oty Ninix/%,. xg ln (eife,)) . 2 dx, 


*#, und *, treten in diesen Gleichungen ganz symmetrisch auf. Man eliminiert 
nun 0f)/0*y, indem man dj,/dx, und df,/dx, auf die andere Seite bringt und die 
erste Gleichung durch die zweite dividiert: 


oy inet oe ig fie: 

CE ee eee (5) 
ata i, eat 1s nid eae 

dx» #o Ho 


Nun beriicksichtigt man, dass die Indizes an den f nur formal sind (f, = fs) 
und zieht die Aussage a) heran, wonach f(¥) noch normierbar ist. Wir diirfen also 
auch die Ableitung zum Beispiel in dem nichtsingularen Punkt 1 vorschreiben. 
Setzen wir also fiir ¥, = 1, /’(1) = A. Dann folgt: 

? Inv + (1/¥) — 1 
f’(*) = A- a, (6) 


= Uline ny 


und daraus durch Quadratur: 


ay i teat 


Dies ist die allgemeine Lésung, wenn man Stetigkeit verlangt. 
Es bleibt noch das Verhalten von f(#) im Unendlichen zu untersuchen, das 
heisst, ob dort das Integral in (7) konvergiert. Hierzu bildet man 
o WE = AE) ee 


é ~€— (mF +4) 
und priift, ob dies fiir § + oo endlich bleibt. Man rechnet leicht nach, dass dies 
nicht der Fall ist; das Integral divergiert also. f’(”) geht etwas schwacher als 1/x 
gegen Null, und ix ) geht fiir x > oo etwas starker gegen Unendlich als Inx. 


yar ae ee | ee 
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Das Integral in (7) kann numerisch ausgewertet werden; durch bereits ver- 
tafelte Funktionen diirfte es sich kaum geschlossen ausdriicken lassen. 

Die Lésungsmethode ist auf andere Mittelbildungen x,y = *%);(%,, ¥,) sowie 
andere lineare Funktionen F(f,) verallgemeinerungsfahig. 


Summary 


The functional equation 


2 HS) = He) + £09) 


occuring in a heat conduction problem is solved analytically. The continuous 


solutions are discussed. 


(Eingegangen: 8. Oktober 1953.) 
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Ludwig Prandtl 


Am 15. August 1953 verstarb in G6ttingen nach langerer Krankheit im 
79. Lebensjahre Lupwic PRANDTL. Sein Name wird unverganglich in die Ge- 
schichte der technischen Wissenschaften eingehen; es gibt fast kein Gebiet der 
angewandten Mechanik, in das PRANDTL nicht fundamentale Fortschritte ge- 
bracht hatte. Von Hause aus Ingenieur, hat er nie verleugnet, dass das wirklich- 
keitsnahe, anschauliche Denken ihm wesentlich naherlag als formale Rechnungen. 
Das ging so weit, dass er sich zum Beispiel nur wenig fiir Minimalprinzipien, wie 
etwa das Hamiltonsche, erwarmen konnte; er wollte die wirkenden Krafte vor 
sich sehen, und wie Lord KELVIN erblickte er in d#/dt nicht einen Differential- 
quotienten, sondern eine Geschwindigkeit. Wenn immer méglich machte er 
kleine Tintenskizzen etwa der zu erwartenden Stromlinien und versuchte dann, 
die Beschleunigungen und Krafte daraus naherungsweise zu entnehmen. Diese 
Skizzen (oft auf alten Briefumschlagen) haben heute bereits eine Art Sam melwert 
erhalten, findet man doch mehrfach darin Anfange wichtiger neuer Entdek- 
kungen. Mit solch minutidser Detailbetrachtung, die immer wieder durch Be- 
ricksichtigung von Versuchsergebnissen korrigiert und erweitert wurde, ging eine 
aussergewOhnliche Kunst im Herausfinden der wesentlichen Ziige parallel. 
PRANDTL verstand es, von den komplizierten Vorgangen gerade so viel wegzu- 
lassen, dass der Rest mit verhaltnismdssig einfachen mathematischen Mitteln 
noch behandelbar war und trotzdem das entscheidende Merkmal erhalten blieb. 

Da ist die beriihmte Grenzschichttheorie, die aus der Beobachtung von Stré- 
mungen in einem sehr bescheidenen Apparat hervorging. (Freilich: dieser Apparat 
enthielt schon die Eckenumlenkschaufeln, die man heute in Windkandlen, Ven- 
tilatoren, Gasturbinen usw. so viel sieht.) Wahrend die gleichzeitige Berticksich- 
tigung von Tragheits- und Reibungskraften im allgemeinen sehr grosse rechne- 
rische Schwierigkeiten bietet, erkannte PRANDTL, dass bei grésseren Reynoldsschen 
Zahlen und guten Formen das Problem aufgespalten werden kann: die Aussen- 
str6mung hat Potentialcharakter, und sie pragt der Reibungs- oder Grenzschicht 
den Druckverlauf auf, so dass deren Verhalten fiir sich allein bestimmt werden 
kann. Damit wares unter anderm auch méglich, die Ablésungen zu erklaren und 
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teilweise vorauszuberechnen. Im Lauf der Zeit entwickelte sich die Grenzschicht- 
theorie ausserordentlich; laminare und turbulente Grenzschichten wurden unter- 
schieden und das Stabilitatsproblem gerade auch von PRANDTL gefOrdert, einer 
befriedigenden Lésung entgegengefiihrt. 

Die Tragfliigeltheorie ist aus einer 4hnlichen Vereinfachung entstanden. Genau 
betrachtet, ist die Umstromung eines endlich breiten Tragfliigels ausserordentlich 
kompliziert. Dadurch, dass man das Problem linearisiert, das heisst das Ge- 
schwindigkeitsfeld der abgehenden Wirbel als kleine Stérung der Hauptbewegung 
betrachtet, gelangt man zu einer linearen Integralgleichung, die heute nach sehr 
praktisch gestalteten Lésungsmethoden in den Flugzeugfabriken routinemassig 
behandelt wird und eine sehr niitzliche Vorausberechnung der Lastverteilungen 
gestattet. Mit dem friih verstorbenen W. BIRNBAUM zusammen hat PRANDTL 
auch den nichtstationaren Fall behandelt und damit die Grundlage fiir die heute 
so wichtige Flatterrechnung geschaffen. Auch gelang die Ubertragung auf Pro- 
peller (BETz). 

Zu den Arbeiten der jiingeren Jahre gehdren Untersuchungen und Berech- 
nungen tiber Gasstv6mungen durch Diisen und in freien Strahlen. Eine bei PRANDTL 
etwas spdter ausgefiihrte Dissertation (STEICHEN) blieb lange Jahre wenig be- 
achtet. Sie enthalt aber die heute so wichtig gewordene Charakteristikentheorie 
der ebenen Uberschallstvémungen, die (zusammen mit BUSEMANN) spater zu einem 
sehr praktischen Rechenverfahren gestaltet worden ist. 

Die Turbulenz hat PRANDTL zeitlebens beschaftigt. Nicht nur wurden in G6t- 
tingen fundamentale experimentelle Untersuchungen tiber Rohrstr6mungen an- 
gestellt (NIKURADSE, SCHILLER), sondern auch wesentliche theoretische Fort- 
schritte erzielt. Es wurde die Austauschlange eingefiihrt und die freie Turbulenz 
behandelt (mit ToLitMieNn und andern), die Stabilitat der laminaren Schichten 
(TIETJENS, TOLLMIEN, SCHLICHTING) untersucht und zum ersten Male eine ratio- 
nelle Erklarung des Umschlages gefunden. 

Schon friih hatte sich PRaNpDTLr fiir die dynamische Meteorologie interessiert. 
Obwohl etwas weniger bekannt, sind seine Arbeiten auf diesem Gebiet von hoher 
Bedeutung geworden. Einer seiner letzten Aufsdtze gibt eine sehr interessante 
Erklarung der heute im Vordergrund stehenden Erscheinungen der Jet-streams. 

Man koénnte noch vom Prandtlschen Staurohr, von den Windkandlen seiner 
Bauart, von Lagertheorie, von Wdrmetibergang in Rohren berichten und hatte 
doch nicht das ganze Wirkungsgebiet PRANDTLsS in der Strémungslehre umfasst. 
Mit vollem Recht hat man ihn als den Vater der neueren Str6mungslehre be- 
zeichnet, doch darf man darob seine weitere Wirksamkeit nicht vergessen. Da 
ist vor allem die Elastizitét und Plastizitét, die ihn oft beschaftigte. Schon seine 
Dissertation iiber Biegeknickung von Balken bezog sich darauf. Bekannt sind 
die Versuche iiber plastische Deformation von Marmor (mit von KARMAN) und 
das Torsions-Seifenhaut-Gleichnis. Die Plastizitat trat in sehr wichtigen Arbeiten 
in den zwanziger Jahren in den Vordergrund (mit NApat). Damals trug er auch 
eine Betrachtung iiber das Kriechen vor, die die Keime des heute so fundamen- 
talen Begriffes der Dislokation enthielt. 

Obwohl er sich grosse Miihe fiir die Vorbereitung seiner Vorlesungen gab, war 
deren unmittelbarer Erfolg, wie er manchmal mit etwas Wehmut zugab, nicht 
sehr gross; PRANDTL fehlte nun einmal das rhetorische Element. Ganz anders in 
Seminarien und bei Diskussionen in kleinem Kreis, da war sein Urteil bis ins 
hohe Alter von unerhoérter Klarheit und Treffsicherheit. In persénlichen und auch 
in politischen Fragen war er riihrend hilflos. Er konnte sich gar nicht vorstellen, 
dass jemand hinterhaltig oder gar bésartig sein konnte. 
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Mit PRANDTL ist ein Gelehrter von saikularer Bedeutung von uns gegangen. 

‘Fast alle seine Arbeiten tragen den Charakter des Definitiven. Man kann sie er- 
_ weitern und ausbauen, aber man braucht sie nicht zu korrigieren. Es diirfte eine 
_ Ehrenpflicht der deutschen Forschung sein, sie in gesammelter Form herauszu- 
_ geben — enthalten sie doch noch zahlreiche Keime zu weiteren Fortschritten. 


J. ACKERE?T 
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: Theory of Elasticity and Plasticity. By H. M. Westercaarp (John Wiley 
& Sons, Inc., New York, 1952). 176 pp., 19 figs. ; $5.00. 


Seit dem Ende des letzten Krieges ist auf dem Gebiete der Mechanik eine 

_ Menge neuer Lehrbiicher erschienen. Vor allem in Amerika hat die Produktion 

an «Texts» eine erstaunliche Breite erreicht, und die zunehmende Tendenz, statt 

wirklicher Lehrbiicher mit klaren Begriffsbestimmungen und prazisen Gedanken- 

_ folgen blosse Aufgabensammlungen von standig wachsendem Umfang zu schrei- 
_ ben, miisste bedenklich stimmen, wenn sie ohne Ausnahmen wire. 

Gliicklicherweise hat uns aber gerade Amerika in den letzten Jahren einige 
der besten Lehrbiicher geschenkt, und dazu darf wohl auch das vorliegende 
Werk gezahlt werden, obschon es leider ein Torso geblieben ist. Der Verfasser, 
der es nach 25jahriger Lehrtatigkeit unternahm, dieses Buch zu schreiben, ist 
vorzeitig verstorben und hat seine letzte Kraft darauf verwendet, wenigstens 

_ die erste Halfte des geplanten Werkes zu vollenden. So ist zwar keine «Theorie 
der Elastizitat und Plastizitat» entstanden, wohl aber eine beispielhafte Einfiih- 
rung in die Theorie der Elastizitat und (von einigen Bemerkungen iiber die Grund- 
lagen der Plastizitatstheorie abgesehen) eine ausgezeichnete Ubersicht iiber die 
wichtigsten Probleme, Methoden und Lésungen der dreidimensionalen Elastizi- 

_ tatstheorie. 

Schon die einleitenden Kapitel, die sich mit dem Ziel und mit der Geschichte 
des behandelten Gebietes befassen, sind in ihrer knappen und prazisen Darstel- 
lung einmalig und wiirden an sich schon die Anschaffung des Biichleins lohnen. 
Die beiden folgenden Kapitel tiber Spannungen, Verzerrungen sowie die Grund- 
gleichungen der Elastizitat und Plastizitat sind weitgehend geometrisch konzi- 

_ piert und kommen damit vor allem dem technisch geschulten Leser entgegen, 
ohne iibrigens auf Strenge zu verzichten. In den letzten zwei Kapiteln wird die 
elastische Feldtheorie entwickelt, wobei die Behandlung mit Hilfe des Spannungs- 
potentials und des Galerkinvektors in den Vordergrund gestellt und an zahl- 

-reichen dreidimensionalen Problemen erlautert wird. Hi. Ziegler 


Vorlesungen iiber theoretische Physik, Band. V: Thermodynamtih und 
Statistik. Von A. SOMMERFELD (Dieterichsche Verlagsbuchhandlung, Wiesbaden 
1952). 374 S., 41 Abb.; geb. DM 24.80, brosch. DM 23.30. 

Der vorliegende Band V, Thermodynamik und Statistik, beschliesst die sechs- 
bandigen Vorlesungen iiber theoretische Physik von A. SoMMERFELD. Es war 
dem Verfasser nicht mehr vergénnt, das Werk selber zu Ende zu fihren; dieser 
letzte Band wurde nach weitgehenden Vorarbeiten SOMMERFELDS von zweien 
seiner ehemaligen Schiiler herausgegeben. Ihnen gebiihrt das Verdienst, die 
Homogenitat des Stoffes gegeniiber den iibrigen Banden und die spezifische Dar- 
stellungsweise des grossen Didaktikers der Physik gewahrt zu haben. 


ZAMP V/12 
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Der thermodynamische Teil ist von betrachtlicher Eleganz infolge der konse- 
quenten Verwendung der «Methode der Potentiale» statt der tiblichen, schwer- 
falligen Kreisprozesse. Viele Anwendungen und Ubungsaufgaben erhellen die 
Niitzlichkeit dieser eigenartigsten Disziplin der Physik, wobei, wie in den 
iibrigen Banden des Werkes, sowohl neueste Entwicklungen (Thermodynamik 
irreversibler Vorgange, um nur ein Beispiel zu nennen) als auch technische Bei- 
spiele zu Worte kommen. 

In der Gaskinetik und statistischen Mechanik wird man vielleicht bedauern, 
dass die Diskussion der Grundlagenprobleme etwas summarisch gehalten ist. 
Dafiir wird man auch hier entschddigt durch eine reiche Auswahl an Problemen 
und Beispielen, die im Text oder als Ubungsaufgaben durchgerechnet werden. 
Als Beispiel sei erwahnt, dass die Berechnung von Warmeleitfahigkeit und Koef- 
fizient der inneren Reibung in Gasen mit Hilfe der Boltzmannschen Integral- 
gleichung bis zum numerischen Resultat durchgefiihrt wird! 

Das Buch reiht sich wiirdig neben die tibrigen Bande und vervollstandigt 
SOMMERFELDS Lehrbuch zum gegenwartig unbestritten besten Lehrbuch der 
klassischen theoretischen Physik, das jeder Physiker besitzen sollte. 

M. R. Schafroth 


Mathematics of Relativity. By G. Y. Rainicu (John Wiley & Sons, New 
York; Chapman & Hall, London, 1950). III + 173 pp.; $3.50. 


Es handelt sich bei diesem Buche um eine in sich geschlossene Darstellung der 
Grundbegriffe der Relativitatstheorie, die von einem Minimum an mathematischen 
Vorkenntnissen ausgeht und den ganzen benétigten mathematischen Apparat in 
vorbildlich einfacher Weise entwickelt; die physikalische Seite der Begriffe und 
Diskussionen wird dabei nie aus dem Auge gelassen, sondern parallel mit der 
mathematischen erértert. — Im einleitenden 1. Kapitel werden im Anschluss an 
die klassische Physik (Mechanik, Elektromagnetismus) die ersten Elemente der 
Tensorrechnung eingefiihrt. Das 2. Kapitel behandelt die Geometrie des vier- 
dimensionalen Minkowskischen Raumes; das 3., daran ankniipfend, die spezielle 
Relativitatstheorie. Im gleichen Verhaltnis stehen die zwei folgenden Kapitel: das 
4. gibt eine sch6ne Einfiihrung in die Riemannsche Geometrie gekriimmter Raume, 
das 5. schliesslich eine Darstellung der allgemeinen Relativitatstheorie, mit den 
einfachsten, aber wichtigen Anwendungen der Gravitationsgleichungen auf be- 
obachtbare Effekte (Planetenbewegung, Kriimmung der Lichtstrahlen, Rot- 
verschiebung der Spektrallinien). Der Verfasser hat auf Vollstandigkeit verzichtet, 
es aber dafiir verstanden, alles Wesentliche streng und in bewundernswerter Klar- 
heit bis zur letzten Einzelheit zu erlautern; das Buch kann als Einfiihrung in 
die Relativitatstheorie warmstens empfohlen werden. B, Eckmann 


Uber Kurven und Flachen in allgemeinen Raumen. Von P. FINSLER 
(Verlag Birkhauser AG., Basel 1951). 160 S.; sFr. 15.40. 


Die nachstehenden Zeilen sind keine «Rezension», sondern ein «Hinweis»: 
das in der Uberschrift genannte Buch ist namlich der unverdnderte (photo- 
mechanische) Abdruck der Géttinger Dissertation von P. FINSLER aus dem 
Jahre 1918. Diese Dissertation bedeutete — neben den gleichzeitig entstandenen, 
durch die allgemeine Relativitatstheorie angeregten grossen Arbeiten von LEVI- 
Crvira, Wry und anderen — einen der wesentlichen Schritte von der «klassi- 
schen» zur «modernen» Differentialgeometrie. Sie kniipfte an Begriffe aus der 
Variationsrechnung an, wie sie besonders von CARATHKODORY in den Vorder- 


: 
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grund gestellt wurden: wahrend in der Riemannschen Geometrie das Langendif- 
ferential einer Kurve durch ds = (2X g,,, dx; dx,)1/2 gegeben ist, wobei die &;, Funk- 
tionen der Ortskoordinaten %,, %5, ..., x, sind und die quadratische Form positiv 
definit ist, setzt FINSLER viel allgemeiner: dsi= F(#,, ...,% 3 044). dx,), wobei 
die Funktion F positiv homogen vom ersten Grade in den dx,, im tibrigen aber (bis 


_ auf Regularitatsannahmen) ganz willkiirlich ist. Auf diesem Grundbegriff wird 


_nun eine dem Vorbild der elementaren Differentialgeometrie folgende Geometrie 
der Kurven und Flachen aufgebaut — iibrigens in eindrucksvoller Einfachheit, frei 


_von raffinierten Formalismen. Man weiss, eine wie grosse und immer wachsende 


Wirkung diese Abhandlung gehabt hat — die Theorie der «Finslerschen Raiume» 
bildet heute ja einen grossen und wichtigen Teil der Geometrie. Die Dissertation, 
die den Grundstein dieser Theorie bildet und auch heute noch lesens- und stu- 
dierenswert ist, war aber nicht in einer Zeitschrift erschienen und nicht im Buch- 


_ handel erhaltlich; darum ist die nunmehr vorliegende Neuausgabe aufs warmste 


zu begriissen. — Neu hinzugefiigt ist ein (von H. SCHUBERT bearbeitetes) sehr 
ausfiihrliches Literaturverzeichnis (schatzungsweise tiber 600 Biicher und Ab- 
handlungen), das sich allerdings nicht nur auf die eigentliche « Finslersche», son- 
dern auf die «allgemeine» Differentialgeometrie im weiteren Sinne bezieht. 
H. Hopf 


Physikalisches Taschenbuch. Von H. EseErt (Verlag Friedr. Vieweg & Sohn, 
Braunschweig 1951). 522 S., 147 Abb.; DM 14.80. 


Wenn man im allgemeinen solch «kleinen Handbiichern» gegentiber eher kri- 
tisch gesinnt ist, bedeutet dieses Buch eine wirkliche Bereicherung aller beste- 
henden Nachschlagewerke, wie sie der Physiker und der Ingenieur im Laborato- 
rium taglich brauchen. In systematischer Folge sind von den Mitarbeitern dieses 
Taschenbuches die Grundgesetze und ihre Anwendungen in verschiedenen Ge- 
bieten der Physik zusammengestellt. Behandelt werden: Mechanik, Akustik, 
Optik, Warme, Elektromagnetismus, Atom und Kernphysik, Werkstoffe und 
die Grundbegriffe der Astrophysik. Besonders hervorzuheben ist das einleitende 
Kapitel iiber die mathematischen Hilfsmittel, enthalt es doch in guter Darstel- 
lung vielgebrauchte Formeln und Satze der Integral- und Differentialrechnung. 

Das «Physikalische Taschenbuch» ist vielleicht gerade deswegen so praktisch 
und niitzlich, weil es neben reinen Zahlenangaben immer einen kurzen und klaren 
Uberblick zu den betreffenden Wissensgebieten vermittelt. O. Huber 


Typical Means. By K. CHANDRASEKHARAN and S. MINAKSHISUNDARAM 
(Oxford University Press 1952). 139 pp.; 35s. 

Seitdem M. Riesz vor mehr als vierzig Jahren die «Typical Means» oder 
«Rieszschen Mittel» zur Summation divergenter Reihen eingefiihrt hatte, haben 
sich diese zu einem kraftigen Hilfsmittel entwickelt. Friihere Darstellungen dieser 
Summationsmethode finden sich in dem Buche von Harpy und Riesz The 
General Theory of Dirichlet’s Series (Cambridge 1915), und im Heft 51 (1931) der 
«Mémorial des Sciences mathématiques» von E. KoGBETLIANTZ. Seither hat sich 
um die «Typical Means» eine ausgedehnte Theorie aufgebaut, zu der namhafte 
indische Mathematiker Wesentliches beigetragen haben, und es ist sehr zu be- 
griissen, dass nun eine Gesamtdarstellung erfolgt ist. Die ersten zwei Kapitel 
behandeln die Summationstheorie divergenter Reihen mit Hilfe der Rieszschen 
Mittel, die zwei letztern Kapitel deren Anwendung auf Dirichletsche und Fourier- 


sche Reihen. A. Pfluger 
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Einfiihrung in die Viskosimetrie und Rheometrie. Von H. UMSTATTER 
(Springer-Verlag, Berlin, Géttingen, Heidelberg 1952). 152 S., 106 Abb.; DM 16.50. 


Die Lehre von der Fliissigkeitsreibung hat sich im Laufe der letzten Jahre, 
angeregt durch die Bediirfnisse der Praxis, zur sogenannten «Rheologie» ent- 
wickelt, die auch Fliissigkeiten umfasst, die nicht dem einfachen Newtonschen 
Reibungsgesetz folgen. 

Im ersten Teil des: Buches beschaftigt sich der Verfasser mit der Theorie 
solcher Korper, die wesentlich auf dem (schon von MAXWELL eingefihrten) 
Begriff der Relaxation beruht. Der elastische K6rper und die Newtonsche Flis- 
sigkeit sind Grenzgesetze fiir unendlich lange bzw. sehr kurze Relaxationsdauer. 

Im zweiten Teil werden die Strémungsformen ausschliesslich fiir zihe New- 
tonsche Fliissigkeiten fiir verschiedene Konturen und Grenzbedingungen behan- 
delt. Fiir Zahigkeitsmessungen sind die kleinen Reynolds-Zahlen natirlich be- 
sonders wichtig. 

Auf dieser Grundlage werden dann im dritten Teil sehr zahlreiche absolute 
Viskosimeter ausfiihrlich beschrieben und wertvolle Winke fiir deren Handhabung 
gegeben. Die Relativinstrumente, die in der Praxis doch haufig verwendet wer- 
den und die natiirlich einer Eichung bediirfen, werden nicht behandelt. Auch ist 
es schade, dass bei den absoluten Instrumenten keinerlei Vergleichsmessungen 
angefiihrt werden. — Dass, wie auf Seite 80 steht, «die Ursache der Entstehung 
der (Prandtlschen) Grenzschichten bisher trotz grésster Bemtihungen nicht auf- 
geklart werden konnte», mutet etwas merkwiirdig an; unter andern hatte das 
eine Zeile vorher zitierte Buch von SCHLICHTING dem Verfasser eine durchaus 
befriedigende Erklarung geben kénnen. J. Ackeret 


Theory of Perfectly Plastic Solids. By W. PraGeErR and P. G. HopcE (John 
Wiley & Sons, New York; Chapman & Hall, London, 1951). 264 pp., 90 figs. ; $5.50. 


In der von I. S. SOKOLNIKOFF herausgegebenen «Applied Mathematics Series» 
smmd"in den-tetzten” fahren eine Reihe sefr bemerkenswerter Monographien aus 
dem Gebiete der angewandten Mathematik erschienen. Dazu gehért auch das 
vorliegende Werk, das aus der in rascher Entwicklung befindlichen Plastizitats- 
theorie ein erstes, heute einigermassen abgerundetes Teilgebiet behandelt. 

Das Buch ist als Einfiihrung konzipiert und elementar geschrieben. Es ver- 
zichtet auf die Exposition der physikalischen Aspekte des Plastizitatsphanomens. 
Es beschrankt sich, wie der Titel sagt, auf perfekt plastische (Mises- und Prandtl- 
Reuss-) Kérper und entwickelt ihre mathematische Theorie unter Betonung der 
technisch wichtigen Fragestellungen, wobei aus theoretischen wie praktischen 
Griinden das ebene Problem im Vordergrund steht. Die einzelnen Abschnitte 
behandeln einfache Fachwerke und Balken, die Torsion von zylindrischen Kér- 
pern, beim ebenen Verzerrungszustand neben axialsymmetrischen und anderen 
speziellen Problemen die allgemeine Theorie samt Grenzanalyse bei beschriankter 
plastischer Deformation, und schliesslich die Extremalprinzipien der Plastizitats- 
theorie. Den Abschnitten sind (teilweise recht anspruchsvolle) Aufgaben bei- 
gefiigt. 

Trotz seines einfiihrenden Charakters stésst das Buch in den letzten Abschnit- 
ten bis zu den aktuellsten Problemen vor, an deren Klarung die Verfasser und ihr 
Mitarbeiterkreis entscheidend beteiligt sind. Es bietet unzahlige Anregungen und 
weist den Weg fiir die kiinftige Ausbildung der Ingenieure auf diesem Gebiet. 


Hi. Ziegler 
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Elektromagnetische Wellen 


Eine Einfiihrung in die Theorie als Grundlage fiir ihre An- 
wendung in der elektrischen Ubertragungstechnik 


von Karl Willy Wagner Dr. phil. et Dr. Ing. h. c. wei- 


land o. Professor an der Technischen Hochschule Berlin 


268 Seiten mit 185 Figuren und 22 Tabellen. — Ganzleinen 
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Das Gebiet der elektromagnetischen Wellen hat bis in die 
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Rohtstr6mung mit schwachem Drall 
Herrn Fr. A. WILLERs zum siebzigsten Geburtstag 


Von Lotuar Corratz, Hamburg!), und HENRY GOrTLER, Freiburg 1. Br) 


Rotationssymmetrische Strémungen zaher Fliissigkeiten mit Drallkompo- 
nente um die Symmetrieachse interessieren aus offensichtlichen praktischen 
Griinden im Hinblick auf strémungstechnische Anwendungen; es wird ihnen 
aber auch ein. grundsatzliches theoretisches Interesse entgegengebracht im 
Rahmen der gegenwartigen Bemiihungen um die Erforschung dreidimensiona- 


_ ler-Strémungsvorgange in zaihen Fliissigkeiten, fiir die sie einen bereits auf- 


schlussreichen, aber noch verhaltnismassig einfachen Sonderfall darstellen. Man 
besitzt bisher nur einige wenige strenge Lésungen der Navier-Stokesschen oder 
auch nur — fiir Stré6mungen an Wanden bei grossen Reynoldsschen Zahlen — 


_ der Grenzschichtgleichungen, welche Strémungen dieser Art darstellen. Dar- 
_ tiber hinaus hat man fiir einige nicht so leicht zugangliche Falle Zuflucht zu 


kaum kontrollierbaren Approximationsansatzen nach dem Muster des Pohl- 
hausen-Verfahrens genommen. 

Bereits die Berechnung der stationaren Stromung mit Drall durch ein ge- 
rades Kreisrohr, die noch relativ einfach tiberschaubar zu sein verspricht, bietet 
bei naherem Zusehen erhebliche Schwierigkeiten, da sich der Berechnung noch 


: grdéssere Komplikationen als im bereits approximativ behandelten Fall der 
_ drallfreien Rohreinlaufstr6mung entgegenstellen. Immerhin zeigt sich, dass im 
Falle hinreichend schwachen Dralls — schwach relativ zur Geschwindigkeits- 


komponente in Richtung der Rohrachse und in diesem Sinne linearisierbar in 


- den Bewegungsgleichungen — in erster Naherung eine Kopplung zwischen der 


Drallkomponente und den anderen Str6mungsabweichungen von der sich asym- 
ptotisch hinreichend weit stromabwarts einstellenden Poiseuille-Strémung nicht 
besteht. Der Fall des schwachen Dralls bietet daher einen ersten Zugang. Auf- 
gefasst als der Fall des Absterbens eines Dralls hinreichend weit stromabwarts 
bildet er den Ausgang fiir die iterative Berechnung eines starkeren Dralls weiter 
stromaufwarts. Auch im Hinblick auf diese weitere Zielsetzung und als Voraus- 
setzung fiir deren Durchfiihrung wird im folgenden eine sorgfaltige Unter- 
suchung des Grenzfalls schwachen Dralls durchgefiihrt. 


1) Forschungsstelle fiir praktische Mathematik der Universitat. . 
2) Mathematisches Institut der Universitat, Abteilung Angewandte Mathematik. 
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Die Berechnung des Dralls fiihrt, wie zunachst gezeigt werden soll, auf ein 
Eigenwertproblem einer gewohnlichen linearen Differentialgleichung zweiter 
Ordnung. Deren Theorie einschliesslich des fiir die Anwendung interessierenden. 
Entwicklungssatzes wird sodann dargelegt. Fiir die interessierenden Eigen- 
werte und Eigenfunktionen werden schliesslich numerische Ergebnisse mitge- 
teilt. Die Genauigkeit der numerischen Ergebnisse ist nicht so sehr auf das 
vorerst hier allein behandelte Problem des schwachen Dralls abgestellt, viel- 
mehr bereits auf die hdheren Anspriiche einer kiinftigen Ausweitung der Unter- 
suchung im oben gekennzeichneten Sinne (mit dem Ziel des Anschlusses an 
eine von einem Anfangsquerschnitt — etwa Rohreinlauf — ausgehende anders- 
artige Entwicklung stromabwarts) zugeschnitten. ° 


1. Grundgleichungen. Grenzfall der Poiseuille-Str6mung 


Es seien x die axiale, 7 die radiale und @ die Winkelkoordinate eines Zylin- 
derkoordinatensystems, y = 0 die Achse und r = 7 die Mantelflache des von 
der zihen Fliissigkeit durchstromten Kreisrohrs. Die Strémungsgeschwindig- 
keit habe die Komponenten wu in x-, v in rv- und w in g-Richtung. Es seien 
schliesslich der Druck, @ die konstante Dichte und v die konstante kinemati- 
sche Zahigkeit der Fliissigkeit. 

Bei Voraussetzung eines achsensymmetrischen Strémungsfeldes sind u, v, w 
und ~ von m unabhangig. Die allgemeinen Navier-Stokesschen Bewegungs- 
gleichungen und die Kontinuitatsgleichung vereinfachen sich dann zu den fol- 
genden vier Grundgleichungen unseres Problems fiir die unbekannten Funk- | 
tionen u(x, 7), v(x, 7), w(x, v) und p(x, 7): 


Witty + Vip = —F2 py tes Hyp + <a (1) 
apy 2 ) ) 
UV_ + VV, — <— = Pe tv lee + Ure + 5 — Fah (2) 
Ete COS { ¥4 Ww, w 
UW, + VW, + - Fie al as v (Orn + Wpy + aa a (3) 
v 
Ug + U,+ ik UR (4) 


Bei undurchlassigem Rohr sind die Haftbedingungen 
“4=v=w=0 fir r=7, (5) 
zu erfiillen. Ferner sind geeignete Vorgaben iiber den Strémungszustand in 


einem Anfangsquerschnitt x = x» und ausreichende Stetigkeits- und Differen- 
zierbarkeitsforderungen zu erfiillen. 
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An Stellen hinreichend weit stromabwarts vom Anfangs- (etwa Einlauf-) 
Querschnitt des Rohrs, durch das die Fliissigkeit unter der Wirkung eines 
Druckgefalles in x-Richtung strémt, werden die radiale Komponente v und die 
Drallkomponente w unter dem Einfluss der inneren Reibung abgeklungen sein, 
und man hat asymptotisch 


v>0, w>0 fir *x>0c0. (6) 


_ Nach (2) und (4) folgt asymptotisch 


fp >0, u,—>90 fir 2% >0e (7) 


und somit 


pb~ p(x), u~u(r) fir x >. (8) 


In nullter Naherung fiir hinreichend grosses x gilt somit nach (1) 


Un + Mg = Px , (19) 


Y ov 


Nach (8) miissen linke und rechte Seite dieser Gleichung konstant sein. Man 


- findet daher durch Integration als nullte asymptotische Naherungslésung fiir 
_ hinreichend weit stromabwarts gelegene Stellen unter Beachtung der Haft- 


bedingung (5) 
2 400U, 
u~ U(r) =Ul1— +), p~ Pi =po- “FE, (9) 


das heisst die bekannte Poiseuille-Strémung, eine exakte Losung der Navier- 
Stokesschen Gleichungen fiir das beiderseits unendlich lange Rohr, fiir den 


_ Fall abweichender Anfangsbedingungen in einem Querschnitt x = x) der asym- 


ptotische Endzustand hinreichend weit stromabwarts von diesem Querschnitt. 


2. Ansatz zur ersten Korrektur des Strémungsverlaufs fiir endliches x 


Um das Abklingen der in einem Anfangsquerschnitt x = %) vorgegebenen 
Abweichungen von der Poiseuille-Strémung in erster asymptotischer Naherung . 
gentigend weit stromabwarts von diesem Querschnitt zu ermitteln, machen 
wir den Stérungsansatz 


u=U(r)+u*, v=v*, w=wrt, p= P(x) +P*. (10) 


Einsetzen von (10) in (1), (2), (3) und (4) und nachfolgende Linearisierung hin- 
sichtlich der Stérungsgréssen unter der Voraussetzung |u*|, |v*|, |w*|< Us 


ZAMP V/7 
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fiihrt auf das System: 


Uut + Ut = — +9 {ut + uf + oe ae (Lay 
ug =— FE + vot, + ot + si ~ +3}, (2.1)% 
Uwr=9 {wx ht a}: (3.1) 

ut tort" =0. (4.1) 


i 


Dieses fiir die erste asymptotische Approximationsordnung geltende lineare 
System zerfallt offensichtlich: Die Gleichung (3.1) fiir w* enthalt nicht die 
jeweiligen u*, v* und p*, und umgekehrt hat w* keinen Einfluss auf u*, v* 
und #* in deren Bestimmungsgleichungen (1.1), (2.1), (4.1). In der ersten 
asymptotischen Naherung ergibt sich also noch sn BEE, Beeinflus- ° 
sung des Dralls einerseits und der Geschwindigkeits- und Druckabweichungen 
von der Poiseuille-Strémung andererseits. (Diese Erscheinung ist offenbar nieht 
auf die hier allein betrachtete Rohrstr6mung beschrankt.) 

Setzt man nur die Kleinheit der Drallkomponente w voraus, so kann man 
die Gleichung (2) noch hinsichtlich w linearisieren. Man erkennt, dass dann (1), 
die so vom Glied — w?/r befreite Gleichung (2) sowie (4) von w frei sind und 
als Bestimmungsgleichungen fiir #, v und # daher wie im drallfreien Fall lauten. 
Ein schwacher Drall hat also an keiner Stelle des Rohrs einen Einfluss auf die 
«Hauptstr6mung» u, v, p. Der schwache Drall (schwach vom Anfangsquer- 
schnitt an) iiberlagert sich einfach der unveranderten Hauptstr6mung. Bei be- 
kannten uw, v und # lasst er sich nachtraglich aus der vollen, in w linearen Glei- 
chung (3) berechnen. Die drallfreie Rohrstr6mung ist von PUNNIS!) approxi- 
mativ berechnet worden unter der Annahme, dass am Einlaufquerschnitt 
u = const und v = 0 tiber den ganzen Querschnitt vorgegeben sind. Fiir diesen 
und fiir jeden anderen Fall der Rohrstr6mung hat die im folgenden durchzu- 
fiihrende Berechnung des Verlaufs eines schwachen Dralls Interesse. Auch 
wenn der Drall erst hinreichend weit stromabwarts vom Anfangsquerschnitt 
hinreichend klein ist, interessiert seine Absterbegeschwindigkeit im Vergleich 
zu dem der Stérungsgréssen u* und v* der drallfreien Hauptstr6mung. 

PUNNIs und vor ihm fiir den entsprechenden ebenen Fall der Kanalstr6- 
mung SCHLICHTING”) haben in ihren Rechnungen eine asymptotische Naherung 
stromabwarts mit einer grenzschichttheoretisch berechneten Naherung vom 
Einlaufquerschnitt her im mittleren Strémungsbereich mit erreichbarer Nahe- 
es SED SETECSA ROSSER Die Berechnung der Drallkomponente w in héherer 


1) B. Punnis, Zur Berechnung der laminaren Einlaufstrémung im Rohr, Dissertation (Gottingen 
ak Nicht im Druck erschienen. 


2) H. Scuxicutine, Laminare Kanaleinlaufstromung, Z. angew. Math. Mech. 14, 368 (1934). 
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asymptotischer Approximationsordnung mit dem Ziele des Anschlusses an eine 
_ Strémung mit stérkerem Drall hinter dem Einlaufquerschnitt, wiirde eine 
Neuberechnung der von der Drallkomponente in zweiter und héherer asym- 

ptotischer Ordnung mitbeeinflussten Gréssen u, v und # erforderlich machen. 
Grundsatzliche Schwierigkeiten wiirden sich dabei nicht ergeben. 

In Ubereinstimmung mit den Untersuchungen von SCHLICHTING und PuNNIS 
_ setzen wir voraus, dass die Reynoldssche Zahl Re = Uy 7,/v hinreichend gross 
ist, um eine weitere Vereinfachung grenzschichttheoretischen Charakters zu 
gestatten: Die Veranderung des Str6mungszustands, hier insbesondere des 
Dralls, soll klein sein gegeniiber der Veranderung in radialer Richtung derart, 
dass man in (3.1) das Ghied w*, gegentiber w% vernachlassigen kann. Der nach- 
folgenden -Berechnung karfin man entnehmen, dass im Rahmen unserer Unter- 
suchung der ersten asymptotischen Naherung fiir einen schwachen Drall diese 
Vernachlassigung bereits fiir mittlere Reynoldssche Zahlen (zum Beispiel 500) 
zulassig sein diirfte. 

Die nachfolgende Untersuchung bezieht sich somit auf die Differential- 
gleichung 
y2 
F 


U, yk : if 
we oe (1 ) w* + ve we Bag w* = 0 (11) 


mit zugeh6rigen Rand- und Anfangsbedingungen. 


3. Zuriickfiihrung auf ein Eigenwertproblem 


Es ist fiir das Folgende zweckmassig, die folgenden dimensionslosen Ver- 
anderlichen einzufiihren: 


v Se (Re aay UW ’e (12) 


ee ne day Vey i 

Damit geht (11) iiber in 
Beis?) ee oat swe — wl (13) 
Neben der Haftbedingung auf s = 1 (Rohrwand) ist Beschranktheit der Funk- 
- tion w* auf s = 0 (Rohrachse) zu fordern, was fiir zweimal stetig differenzier- 
bares w* nach (13) mit Verschwinden von w* auf der Rohrachse gleichbedeutend 


ist. (Aus physikalischen Griinden wird man dies von vorneherein fordern.) 
Die Randbedingungen lauten daher 


a) furs =O und fir s= 1. (14) 
Hinzu kommt die Vorgabe von w* in einem Anfangsquerschnitt, der ohne Ein- 
schrankung mit g = 0 bezeichnet werden kann: 


w*=h(s) fir g=0. (15) 
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Das Verschwinden von w® fiir g oo (Absterben des Dralls unter dem Einfluss 

der inneren Reibung) muss jeder Lésung von (13) und (14) bei beliebiger An- 

fangsverteilung (15) als Eigenschaft zukommen. 
Der Separationsansatz 


w* = F(q) G(s) , (16) 
fiihrt bei Einsetzen in (13) auf 
F(q) = e*# (17) 
und liefert fiir G(s) die Differentialgleichung 
s?G"+sG'—G=i/3s? (rs) ee (18) 


wobei Striche Ableitungen nach s anzeigen und 4 Separationsparameter ist. 
Die Randbedingungen (14) gehen iiber in 


Gi =.G(1} oe0. (19) 


Es wird im nachsten Abschnitt nachgewiesen werden, dass das Eigenwert- 
problem (18), (19) abzahlbar unendlich viele Eigenwerte i, und zugehdérige 
Eigenfunktionen G,(s) besitzt, dass ferner die Eigenwerte samtlich negativ sind 
(in Ubereinstimmung mit dem zu erwartenden Absterben des Dralls mit wach- 
sendem g) und dass unter gewissen Voraussetzungen die Anfangsverteilung 
(15) nach den Eigenfunktionen entwickelbar ist: 


h(s) = 3c, G,(s) . (20) 


v=1 


Die gesuchte Drallverteilung in 0 <s <1, 0 Sq ist dann 


= J c, e7#4 G(s) (21) 


vel 


oder in den urspriinglichen Variablen 


w*(x, 7) -s c, eXp (4) G+). (21a) 


0 


4, Theorie der Eigenwertaufgabe und Entwicklungssatz 
Die in (18) und (19) aufgetretene Eigenwertaufgabe hat die Form 


s?G"(s) + s G'(s) — G(s) = As? f(s) G(s); G(0) = G(1) = 0 (22) 
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mit f(s) =.1 — s*, wobei Striche Ableitungen nach s bedeuten. Die Differential- 
gleichung werde in der selbstadjungierten Gestalt geschrieben: 


L[G]=[s G(s)’ — — G(s) =A f(s) G(s). (23) 


Die allgemeine Eigenwerttheorie bei Differentialgleichungen ldsst sich wegen 
der Singularitaéten bei s = 0 und s = 1 [bei s = 1 verschwindet /(s) = 1 — s?] 
nicht unmittelbar anwenden; aber man kann der Aufgabe eine dquivalente 
Integralgleichung zuordnen. Die Randwertaufgabe 


L[G]= 1s); G(0) = G(1) =0 (24) 


besitzt die mit Hilfe des Fundamentalsystems s, 1/s von L[G] = 0 leicht auf- 
stellbare Greensche Funktion H(s, #): 


yds 


if 
a7 8 bie (ORS Sy 
TTS: Bucs vey 
tf bey Gis, See 

ZS! 


welche die Randwertaufgabe (24) durch 
| ae : 
Bins J H(s, t) r(t) dé (25) 
0 
lést; somit geht (22) uber in 
1 
Goes i t f(t) H(s, t) G(t) de. 
0 


Im folgenden setzen wir /(s) als in <0, 1> stetige, in 0 Ss <1 positive Funk- 
tion voraus. Die speziellen Gegebenheiten der Funktion 1 — s? beriicksichtigen 
wir, indem wir /(1) = 0, //(1) + 0 verwenden. Dann kann man den Kern in 
bekannter Weise symmetrisieren. Mit 


K(s,#) =Vsti(s) ) H(s,); @(s) =Vs f(s) G(s) (26) 
erhalten wir Ss 


os) =A af K(s, 2) y(t) dt. 


Der Kern K(s, ¢) ist stetig und symmetrisch und gestattet die Anwendung der 
Integralgleichungstheorie. Es sei /(s) eine zweimal stetig differenzierbare, bei 
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s=0 und s=1 dreimal differenzierbare, den Randbedingungen h(0) = A(1) 
= k(1) = 0 geniigende Funktion, wobei 


= L{h(s)] (27) 
gesetzt ist. Aus (25) und (26) folgt 


(s) Vs f(s) = | x K(s, 8) ae dt. (28) 


Wenn die Funktion 


in <0, 1> quadratisch integrabel ist, ist#) 
= h(s) Vs {(s) 


quellenmiassig darstellbar, stetig und nach normierten Eigenfunktionen 9,/(s) 
des Kernes in eine absolut und in <0, 1> gleichmassig konvergente Reihe 


ts 


Baise g,(s) mit c,= feo) y,(s) ds (29) 


entwickelbar. Es ist also zu priifen, ob &*(s) in <0, 1> quadratisch integrabel 
ist. Die Stelle s = 0 ist harmlos; es ist nach Definition 

R(s) =sh"(s) + hi(s) — 2S, 
dabei ist h(s) zweimal stetig differenzierbar und /(0) = 0; mithin ist h’(s) — h(s)/s 
in <0, 1> stetig differenzierbar und hat den Grenzwert 0 fiir s > 0; oo hat 
auch k(s) den Grenzwert 0 fiir s > 0 und ist bei s = 0 noch differenzierbar. 
Daher ist k*? an der Stelle s = 0 rechtsseitig stetig (und hat dort den Wert 0). 
Wegen der Voraussetzungen ist k(s) bei s = 1 linksseitig differenzierbar und 
hat dort den Wert 0, so dass k*? auch bei s = 1 (linksseitig) stetig ist; k*(s) 
ist mithin quadratisch integrabel. Aus (29) folgt unmittelbar 


= SCG) anit [s (s) G,(s) h(s) ds , (30) 
wo G,(s) durch 
fs i(s) G2(s) ds = 1 (31) 


1) Vergleiche etwa G. HaMet, Tntegralgleichungen (Springer, Berlin 1937), S. 72, 76. 
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- normierte Eigenfunktionen von (22) sind. Aus der Willkiirlichkeit von h(s) folgt, 

dass der Kern K(s, ¢) nicht entartet ist und unendlich viele reelle Eigenwerte 
_ besitzt. Die Eigenwerte sind samtlich negativ. Schreibt man nadmlich (23) fiir 
eine zum Eigenwert A, gehdérige Eigenfunktion G, an, so folgt durch Multi- 
_ plikation mit G, und Integration iiber das Intervall von 0 bis 1 


J (sors 2 oi) a 
A, = — : as (32) 


v 1 , 


[s#Geds 


0 


_man kann hier G, als reell annehmen und liest 2,< 0 unmittelbar aus (32) ab. 

Zusammenfassung: Die Eigenwertaufgabe (22), wobei f(s) in: <0, 1) stetig, 

in 0 <s < 1 positiv und /(1)-= 0 ist, /’(1) existiert und + 0 ist, besitzt eine 

abzahlbare Folge von saémtlich negativen Eigenwerten 4, mit (reell annehm- 

baren) Eigenfunktionen G,(s). Jede in <0, 1> zweimal stetig und in s = 0 und 
s = 1 dreimal stetig differenzierbare Funktion h(s) mit 


h(0) = (1) = Lfh(s)],.. = 0 


ist in <O, 1> in eine absolut konvergente Reihe (30) nach den Eigenfunktionen 
_ G, entwickelbar. 


5. Numerische Ergebnisse!) 


Die hier gesammelten Erfahrungen werden etwas ausfiihrlicher mitgeteilt, 

als es fiir die vorliegende Arbeit allein n6tig ware, weil die Erprobung der 
-Methoden auch fiir andere dhnliche Faille niitzlich sein kann. 

Die Formen der ersten fiinf Eigenfunktionen in Figur 1 wurden aus den 
Naherungswerten des gewohnlichen und eines verbesserten Differenzenverfah- 
rens mit den Schrittweiten A = 0,2, h = 0,1 und A = 0,05 gewonnen. Die Durch- 
fiihrung auch bei kleiner Maschenweite gestaltet sich sehr bequem und iiber- 
sichtlich, wenn man als Normierung G(1 — h) = 1 benutzt, also von G(1) = 0, 

G(1 — h) = 1 ausgehend, nach den Differenzengleichungen rekursiv G(1 — 2h), 
G(1 — 3h), ... berechnet, wobei man einen geschatzten Naherungswert A fir 
den betreffenden Eigenwert verwendet. Der Wert G(0) wird dann nicht gleich 0 
“sein, sonderner-hangt von A ab. Durch kleine Variation von A und Interpola- 
tion ist es leicht, A so zu bestimmen, dass G(0) = 0 wird. Das Differenzenver- 
fahren ist fiir Berechnung der Gestalt der G,(s) gut geeignet, wahrend zur 


1) Betreffs einer genaueren Darstellung der Methoden vergleiche etwa L. CoLLatz, Eiigenwert- 
aufgaben mit technischen Anwendungen (Akademische Verlagsgesellschaft,” Leipzig 1949). Hier 
“werden im wesentlichen nur die Ergebnisse mitgeteilt. Fiir numerische Hilfs- und Kontroll- 
rechnungen danken wir den Herren J. GRUNscu und M. SiMosseEc. 
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genaueren: Ermittlung der Eigenwerte 2, und zur Aufstellung von Schranken 
fiir die A, andere Verfahren giinstiger erscheinen. 2 


(A 42-21, 38230 ) 


Ommn02 <a 04s CANO caotese 


Figur 1 
Verlauf der ersten fiinf Eigenfunktionen. 


Obere.Schranken fiir die Betrage der A, liefert das Ritzsche Verfahren. Der 
Naherungsansatz 


mit der Forderung 


ie 


=0 mit J/= 7c (s — s8) w2(s) + s w’%(s) + - w°(s) ds (p = lose 


wurde fiir # = 1, 2 und 3 durchgefiihrt und lieferte die oberen Schranken 


Vol. V, 1954 Rohrstrémung mit schwachem Drall 105 


Obere Schranken fiir 


oe) 21,788 154,2 - | 
p=3 21,3986 86,33 587 


An X 2 


Ts eee ; 216 
noo 7 L z : 
/yi-e®a 
0 


(33) 


welche die Naherungen A, ergibt: 
ee ee es | 3 | 4 | 5 


| a 16 | 64 | 144 | 256 | 400 


Das Verfahren der schrittweisen Naherungen liefert bequem obere und untere 
Schranken fiir 2,. Ausgehend von einer stetigen Funktion /)(s) mit 


werden weitere F, nach 
LiF l= (so 8°) 3; -2(0) = 71) = 0 


und Schwarzsche Konstanten 


1 
ay= | (8 = 9) FFs ds) 
0 


’ Quotienten 
ay), 
Pe = | (ha Op 1,03} 
berechnet und die Abschatzungsformel 
a a YE) Ay S Pa 
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benutzt, wobei —/, eine untere Schranke fiir —/, bedeutet mit 


= Psy < ly < Ag - 


Fir 
F(s) = 23+ Ss* 4 os Zs? 
wird 
873470 18296768 _ hy 
He = — 4ogaq = 72138709, bs = mae 21 ,38340 
und damit 


—A, = 21,38266 + 0,00074 (Fehler also héchstens 0,0035%). 


Fir die niedrigeren Eigenwerte ist hier der Potenzreihenansatz 


G(s) = 7a, 8"? (34) 
v=0 
sehr geeignet ; aus der Rekursionsformel 
A 
sae eae (2,2 — My_4) 
lassen sich, von % = 1 (und formal «_, = 0) ausgehend, «%, «4, ... bequem 
berechnen ; man schatzt A, bildet dann 
a 
Ga) = Da, 
vp=0 


und variiert A so, dass ®(A) = 0 wird. Auf diese Weise werden 


—A, = 21,382302, —A, = 74,7793, -—A, = 160,15 
ermittelt. 


Untere Schranken auch fiir einige héhere Eigenwerte kann man aus den 
reziproken Potenzsummen der Eigenwerte gewinnen. Es ist 


ee 1 


ee i: Es [ Ko, s) ds , 


0 


ms 1 Lean 

. -s 1 ames : _ ~ a i f i a y. oO 

j= a 2 a K,(s, s) ds = ie [K(s, t)]2 ds dé, 
( 


ti 
= 3 ia =| / [K,(s, t)]? ds dt . 
: 0 0 
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 Dabei ist der iterierte Kern 
i 
K,(s, 1) = [| K(s, w) K(s, t) du 
é 


benutzt; es ist hier 


K,(s,)=4Vie—s9 ¢- 9) { (3 — 104 5s? — st) — stins } 


fir OStss, K,(t,s)=K,(s,#). 


_ Man erhailt 
eee S020. 089' 3394s 
Ss eee =? 0,002 430555)5,. 1... . 
S=  ogaoagg = 0,000 004 81760. .... 


Bei den in Figur 1 angeschriebenen Naherungswerten fiir die Eigenwerte A, 
und A,;, welche aus sj, Sy, s, in Verbindung mit der asymptotischen Formel (33) 
gewonnen wurden, ist die letzte hingeschriebene Ziffer als unsicher anzusehen. 


6. Anwendung 


Der vom Betrage kleinste Eigenwert 1, = — 21,38230 liefert das in Hauptstr6- 
mungsrichtung am langsamsten abklingende und daher stromabwarts schliess- 
lich dominierende Glied in der Reihe (21a) fiir die Drallkomponente w*. (Dies 
gilt fiir den im allgemeinen bestehenden Fall c, + 0.) Aus der Untersuchung 
von Punnis!) ist zu entnehmen, dass: in den entsprechenden Entwicklungen 
fiir u* und v* der vom Betrage kleinste Eigenwert zu — 16,03966 berechnet 
wurde (dort ohne Angabe des Genauigkeitsanspruchs), und dass daher das in 
diesen Entwicklungen hinreichend weit stromabwarts schliesslich tiberwiegende 
' Glied wie exp (— 16,03966 x/7) Re) abklingt. Wahrend dieses Glied also auf 
- einer Strecke der Lange 7, Re/16,03966 auf das (1/e)-fache abnimmt, geschieht 
dies bei dem in der Entwicklung von w* schliesslich dominierende Glied auf 
einer Strecke der Lange 7) Re/21,38230, somit bereits nach Durchlaufen von 
drei Vierteln der erstgenannten Strecke. 

Um nun an einem Beispiel etwas eingehender das Absterben eines (schwa- 
chen) Dralls in der Rohrstrémung zu verfolgen, seien im folgenden numerische 
Ergebnisse der Anwendung auf den Fall mitgeteilt, dass im Anfangsquerschnitt 


1) B. Punnis, Zur Berechnung der laminaren Einlaufstromung im Rohr, Dissertation (Gottingen 
1947). Nicht im Druck erschienen. 
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qg=0 die Drallverteilung w*(0, s) = h(s) =es (€ = const, |e|< Up) vorge- 
geben ist. Eine Realisierung dieses Falls kann man sich so denken, dass die 
Fliissigkeit vor Eintritt in den Querschnitt g = 0 ein hinreichend langes, mit 
der konstanten Winkelgeschwindigkeit w = e¢/r, rotierendes Rohr durchlaufen 


10 


Figur 2 & 
Verteilung der dimensionslosen Drallgeschwindigkeit tber den Rohrradius fiir verschiedene Quer- 
schnitte g = const. 


hat, so dass in diesem Rohr vor seinem Endquerschnitt ¢g = 0 schliesslich eine 
Strémung besteht, die aus einer Uberlagerung der Poiseuille-Stro6mung in 
Achsenrichtung und der starren Rotation zusammengesetzt ist. Diese Str6mung 
geht mit w*(0, s) = es im Querschnitt ¢g = 0 in unser ruhendes Verlangerungs- 
rohr von gleichem Durchmesser iiber. 

Der Tatsache, dass die in Abschnitt 4 gestellte Forderung h(1) = 0 durch 
die angenommene Anfangsverteilung nicht erfiillt ist, kann man durch eine 


A 
; 


4 
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-numerisch nicht ins Gewicht fallende Modifikation Rechnung getragen denken: 
Man gebe h(s) = es in 0 Ss <1 — 6 mit einem beliebig klein zu denkenden 
-positiven 6 vor und setze h(s) stetig und zweimal stetig differenzierbar so fort, 
dass es zwischen 1 — 6 und 1 steil abfallend auf 4(1) = 0 abnimmt. (Praktisch 
-gedeutet begeben wir uns aus dem Quersohnitt des unstetigen Ubergangs be- 

liebig wenig stromabwarts, wo sich wegen des Haftens an der ruhenden Rohr- 
_wandung bereits eine noch sehr diinne Grenzschicht der Drallkomponente aus- 
-gebildet hat.) Diese mit Riicksicht auf die Theorie erforderliche Modifikation 


Figur 3 
'-Verteilung der Drallgeschwindigkeit im Rohr mit Ort des Maximums in den einzelnen Querschnitten. 


hat bei hinreichend kleinem 6 keinen Einfluss im Rahmen der Genauigkeit der 
numerischen Auswertung und braucht bei dieser auch nicht beriicksichtigt zu 
werden. oj 

Die Auswertung erfolgte unter Beriicksichtigung der aus Abschnitt 5 nume- 
risch bekannten fiinf ersten Eigenwerte und Eigenfunktionen. Wahrend diese 
zu einer guten Approximation der Anfangsverteilung noch nicht ausreichen, 
zeigt sich, dass bereits in g = 0,01 das fiinfte Glied der Reihe (21a) nur noch 
maximal einen Beitrag von 4% des Ortlichen w* liefert [in etwa s = 0,1, vgl. 
G;(s) in Figur 1]. Im Querschnitt g = 0,06 erreicht sogar das zweite Glied der 
Reihe nur noch einen Beitrag von maximal etwa 5% von w*, und weiter strom- 
abwarts ist praktisch nur noch die erste Eigenfunktion G,(s) fiir die Gestalt 
der Drallverteilung tiber den Rohrquerschnitt massgebend. Von hier ab sind 
also die Drallverteilungen in den Querschnitten ahnlich, und fiir die Abnahme 
des Drallbetrags ist nur noch der erste Eigenwert A, bestimmend. Insbesondere 
reicht fortan die wandnahe Zone, in der der Drall w* von seinem Maximum 
im betreffenden Querschnitt monoton auf Null abnimmt, von etwa s = 0,43 
[Ort des Maximums von G;,(s)] bis s = 1. 
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Figur 2 zeigt die Verteilung des Dralls, genauer: w*/e, tiber den Rohr- 
radius fiir eine Folge von Querschnitten g = const (q = x »/76 U,). Das Her- 
einwandern des Drallmaximums von der Rohrwand bis zur asymptotischen 
Endlage (etwa s = 0,43, siehe oben) ist in Figur 3 nochmals veranschaulicht, 
zusammen mit den an ihren Querschnittsorten aufgetragenen Drallverteilungen. 


400 


0 0,04 0,08 012 016 020 
x:Vv 
Ae 
Figur 4 I U, 


Verteilung der Winkelgeschwindigkeit des Dralls auf der Rohrachse. 


Von der Rohrachse s = 0 aus wachst der Drall in den einzelnen Quer- 
schnitten zunachst von Null praktisch linear an. Grenzt man den achsennahen 
Bereich ab, in welchem sich der Drall in den einzelnen Querschnitten mit vor- 
gegebener Genauigkeit durch die Tangente in s = 0 approximieren lasst, so 
erhalt man eine Kurve, die analog zur Maximumskurve der Figur 3 verlauft, 
raumlich also ein sich in Strémungsrichtung verjiingendes Rohr. In jedem 
Querschnitt hat der Drall innerhalb dieses Rohres approximativ die Verteilung 
einer starren Rotation; von Querschnitt zu Querschnitt stromabwéarts klingt 
deren Winkelgeschwindigkeit [0w*/0s],_» jedoch rasch ab (vgl. Figur 4). 


Summary 


The steady laminar flow at high Reynolds’ numbers in a straight tube of cir- 
cular cross-section is investigated under the condition that swirl is present. Given 
any axially symmetrical distribution of the rotational velocities, the flow down- 
stream is calculated in the special case that the rotational velocities are small 
compared to the axial velocity. The problem is reduced to an eigen-value 
problem of an ordinary linear differential equation of second order, the theory 
of which is developed in detail. Numerical values for the five lowest eigen-values 
and the corresponding eigen-functions are given. As an example, the decay of 
an initially rigid rotation is considered. 

(Eingegangen: 16. Juni 1953.) 
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The Role of Coordinate Systems in Boundary-Layer Theory’) 


By Sau Kapiun, Pasadena, California?) 


1. Introduction 


For a given problem, the approximate solution given by boundary-layer 


_ theory depends on the system of coordinates used when the simplifying assump- 
_ tions of boundary-layer theory are applied to the Navier-Stokes equations. In 


general, different systems of coordinates lead to boundary-layer equations which 
are not equivalent, that is, their solutions represent different flow fields. The 


; object of this paper is to investigate in detail how these flow fields depend on 


the choice of coordinates, how different boundary-layer approximations to the 


_ same flow problem are related to each other, to the external flow and to the 
_ flow due to displacement thickness. The discussion is restricted to incompres- 
_ sible, steady, two-dimensional flow without separation. However, many of the 
_ results hold much more generally and will be discussed in a later paper. 


The main result of this paper is contained in Theorem 2. Normally, one uses 
boundary-layer theory in the following way in order to obtain a picture of the 


- complete flow field of a viscous fluid (outside the wake) : The flow field is divided 


into two separate regions, that is, a boundary-layer region where the flow field 
is obtained from boundary-layer equations, and an outer region where the 
Euler equations are used to obtain an external flow, corrected for the displace- 
ment effect of the boundary layer. There has been considerable discussion about 
where and how to patch the two parts of the flow field, and about how to 
proceed to higher order approximations. However, Theorem 2 shows that one 
can find a system of coordinates, such that the boundary-layer solution with 
respect to this system gives an approximation which is valid in the whole flow 
field. Both the external flow and the flow due to displacement thickness are 
included analytically in this approximation and, hence, the problem of patching 
is automatically eliminated. A coordinate system with these properties will 
be here referred to as optimal. In Section 6, the problem is discussed to what 
order such an optimal boundary-layer solution is valid. In general it gives a 
better approximation to the exact flow field than does the composite flow field 
described above. It also appears to form a reasonable starting point for finding 


1) Research prepared under Office of Naval Research Contract N 6 ONR-244, Task Order VIII, 


Project NR 061—036. 
2) Guggenheim Aeronautical Laboratory, California Institute of Technology. 
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higher order approximations, that is, approximations which are definitely out- 
side the scope of boundary-layer theory. 

In order to find such an optimal system of coordinates it is in general first 
necessary to compute the external flow and the flow due to displacement 
thickness. But once this system has been found, the optimal solution can be 
found immediately by means of another theorem (Theorem 1), which states 
that the boundary-layer solution with respect to one coordinate system, say an 
optimal one, may be found by a simple substitution into the solution with 
respect to any other system, say a conventional one. 

Concrete examples illustrating the application of the above general theorems 
are given in Section 5. 

In this paper, the concepts of the boundary-layer solution, and of the solu- 
tions for external flow and for flow due to displacement thickness, are defined 
with the aid of two special limiting processes. This viewpoint has previously 
been adopted by several authors, for example by H.WEyYL[1]1), K.O. FRIEDRICHS 
[2], and G. E. Latta [3]. This method is actually only a formal and more precise 
restatement of PRANDTL’s original method, and hence is not widely used. 
However, it appears very natural for deriving the results of the present paper 
and the reasoning will be based on systematic use of the two limiting processes. 

There are many examples in the literature of the use of various systems 
of coordinates in connection with boundary-layer theory and, in particular, it 
has been noted that they give the same skin friction and agree approximately 
within the boundary-layer region proper. To the author’s knowledge, however, 
no systematic study of the relation between the solution based on the various 
systems has been undertaken and, especially, the two theorems referred to 
above appear to be new. The starting point for the present investigation was 
actually the following result due to M. D. VAN Dyke [4]. Consider Oseen flow 
past a semi-infinite flat plate. If rectangular coordinates are used to make a 
boundary-layer approximation, the resulting approximation resembles essen- 
tially the Blasius solution, and is valid only in the boundary layer proper. On 
the other hand, if parabolic coordinates are used, the boundary-layer approxi- 
mation satisfies exactly the full Oseen equations (cf. Example 2 in Section 5). 
This result was further analyzed by P. A. LAGERSTROM, who pointed out how 
in other special cases coordinates may be found which are optimal in the sense 
described above [5]. He also suggested to the author the general problem of 
finding optimal coordinates, and guided the present study. The ideas developed 
by G. E. Latta [3] have also influenced the author. A closer comparison of [3] 
with the results and methods of the present paper would undoubtedly shed 
further light on the problems studied here. This comparison has, however, as 
yet not been carried out. 

A preliminary announcement of the results of this paper was given in [6]. 


1) Numbers in square brackets refer to References on page 134, 
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2. Boundary-Layer Approximations as Limits of Exact Solutions 


The Navier-Stokes equations. The equations for steady, viscous, incompres- 
sible flow are 


(q-V)4=-Vet+rp'q, (1a) 
V-q=0, (1b) 
where » is the kinematic viscosity and # is the pressure divided by density. 


In two-dimensions, the kinematical form of (1) is obtained by taking the 
_ curl of (1a) and integrating (1b) by means of the stream function yp, which gives 


q:‘Vo=r7f7a, (2a) 
o=—p*»y, (2b) 

0 0 
es ya ——t, (2c) 


where w is the vorticity. 
In arbitrary curvilinear coordinates, &*, (€1 = &, €? = n) the corresponding 
equations can be expressed as 


wi DH = —|g| gt 2% 4 oV[e| gt eh (1a’) 
ee = 0, (1b’) 
and 
Coma ee ge, 2 
o) Tia _ (Viel ght oe), (2b’) 
ua 5, p= 5h. (2c’) 


Here y, @, and v have been treated as absolute scalars. The velocity compo- 
nents u and v may be considered defined by (2c’). In the language of tensor 
analysis [7], they are actually the contravariant vector-density components of 
the velocity?). It should be noted that the square of the magnitude of the 
velocity vector is u’ u/ g;;/|g|, so that the value of u* is not equal to the value 
of the projection of the velocity vector on the é*-axis (cf. Section 3 below). 

1) idl sete . es here denote first and second covariant derivatives of the vector-density 

De’ DEE DE 
ut (cf. [7], p. 82). Sees 

2) Wevt [1] restricts himself to conformal coordinates, in which case w and v are also equal 

to the covariant vector components of the velocity. 


ZAMP V/8 
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In this paper we shall be concerned with viscous flow past (or towards) a 
rigid body. The curve representing the surface of the body (called the wall) 
need not be closed, but may, for example, correspond to a semi-infinite flat 
plate or a wedge. It will be required, however, that no other boundaries be 
present. 

The boundary conditions at the surface of the body are then 


t=O 0 (3a) 


at the wall. y must then be a constant at the wall, which will be normalized to 
zero. 

At infinity the customary boundary conditions are (in suitable Cartesian 
coordinates) 


q>Ui, p>, as *>-oo or |y|>oo. (3b) 


However, for sufficiently divergent infinite bodies there is no potential solution 
of the equations for non-viscous flow which satisfies (3b). In this case we 
require that 
eet (3c) 
Yop 
at infinity, where wy, represents one of the solutions for potential flow past the 
body. Evidently (3b) is a special case of (3c). 

The term exact solution for flow past the body will henceforth be used to 
denote a solution of (1) or (2) with the boundary condition (3a) satisfied at the 
body and (3b) or (3c) at infinity. More precisely, the term will denote a class of 
solutions with v as a parameter (0 < » < v9) where all solutions satisfy identical 
boundary conditions. Boundary layer theory is concerned with approximations 
to the exact flow field for small values of y. It will be introduced below with the 
aid of two limit processes. 

First limit process. Let f(P, v) be a function of position and viscosity. Here 
P is used as an abbreviation for the coordinates (&, 7) of a point. In the first 
limit process, denoted by lim,, P is held fixed while » tends to zero: The 
resulting limiting function is designated by the subscript e. 

f.(P) = lim, f(P, ») Sinus ”). (4) 
If this limit process is applied to the exact solution for flow past a solid, that is, 
to y(P, v), w(P, v), and p(P, v), the resulting limiting functions represent a 
flow which will be referred to as the external flow. 

In applying the first limit process to equations (1) and (2) at a point where 
the derivatives appearing in (1) or (2) actually tend to the corresponding deri- 
vatives of the limiting functions (which may be called a regular point), one may 
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simply replace py by y,, etc., and put » = 0 in the equations. The resulting 
equations are the Euler equations, i.e., the equations of motion of a perfect 
fluid. 


Applying the first limit process to the boundary conditions, however, it is 


- seen that, since uw = v = 0 at the wall, w, and v, must be zero at the wall, while 


A 
5 
a 


Y/Y» > 1 at infinity. There is, in general, no solution of EULER’s equations, 


satisfying the Euler equations at every point of the flow, which also satisfies 
these boundary conditions, so that all points of the flow cannot be regular in 
the above sense. 

While many types of irregular behavior are conceivable, the type made 


_- prominent by experience consists of a discontinuity in the tangential velocity 


(a vortex-sheet discontinuity), which represents a line of slip in the perfect 
fluid. The behavior under the first limit process at such a discontinuity is not 
regular. As »y + 0 at a point on the line of discontinuity, the viscous terms in 
(1) or (2) do not tend to zero, and EULER’s equations are not approached in the 
limit. Once such lines of slip in the perfect fluid are permitted, EULER’s differ- 
ential equations are not everywhere satisfied, and an infinite number of 
perfect-fluid flows becomes possible. The relevant one, that is, the limit which is 
actually approached, is decided by the action of viscosity, in a manner which, 
as yet, is not understood in many essential respects. 

In ordinary problems, experience and physical arguments show that the line 
of slip occurs right at the wall, following the wall for some distance downstream 
from the point of impingement of the streamline coming from upstream infinity, 
but that it may eventually separate from the wall, defining a wake. It is to be 
expected that the limits y,(P), ui(P), p,(P), actually exist and satisfy EULER’s 
equations, at least outside the region of the wake, which is essentially unex- 
plored at the present time?). The present discussion will, however, be restricted 
to flows where no separation occurs. Then, all the streamlines come from up- 
stream infinity, the flow is irrotational and satisfies LAPLACE’s equation. The 
action of viscosity is restricted to the determination of circulation about the 
body, in accordance with the Kutta condition. y, can then be found by solving 
the corresponding potential problem. 

If the relevant solution of EULER’s equations, uj, is used as a starting point, 
the only essential effect of introducing y small, but greater than zero, is to 
replace the discontinuity by a rapid, but continuous transition, which is nearly 
completed within a narrow region (boundary layer). Perturbation problems 
of this nature are referred to in recent literature as singular (see [3] and referen- 
ces given there). 

In a regular perturbation problem, the successive perturbations are obtained 
by successive application of the first limit process, which is essentially equiva- 
lent to expanding the solution in powers of a small parameter. In a singular 


1) This is briefly commented on further in Section 6. 
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perturbation problem, however, u*—> ui non-uniformly in any region contain- 
ing boundary points. Hence, at » >0, the limit uj breaks down as an approxi- 
mation, sufficiently close to the boundary. An approximation for the con- 
tinuous transition across the boundary layer may be found by the application 
of another limit process, described in the next paragraph. 

Second limit process. Let ¢ be a small parameter which, in a sense, measures 
the effective thickness of the transition region. In the present case e will be 


taken proportional to Vv. More generally, » may be a function of ¢ such that 


lim 7 = const + 0. (4b) 


e—0 


A system of coordinates (€, 7) will be chosen with the only restriction that 7 
be zero at the wall. It is customary also to require (cf., for example, [8], p. 119) 
that the lines € = const be orthogonal to the wall. This assumption will not 
be made here. The wider choice of coordinate systems thus provided will be 
essential for obtaining the principal result of this paper. 

In the second limit process, the point P will not be fixed but will move 
toward the wall along a curve € = const in such a way that the ratio 


n= (5) 
is fixed, in other words, keeping € and 7 const. Let f(P, v), as before, be a 
physical quantity depending on position and viscosity. The limiting function 
approached in the second limit process will be denoted by f;(€, 7) or lim; f and 
defined by 

f(E, ) = lim. f= lim 7S. n, é) . (6) 


e—>0;6,7 fixed 


Here ¢ is used as an abbreviation for the coordinate system (é, 7). Note that 
the first limit of / depends on f only; the second limit, on the other hand, 
depends also on the choice of coordinates. To emphasize this, ¢ is used as a 
subscript in the notation introduced by (6). 
Now consider an exact solution as described by y, uv, v, w, p. We define 
ae geld . 


py ND ee (7) 


The boundary-layer approximations of the quantities , u,v, @, P may then be 
defined as the limiting functions w;, ue, Ue, @e, and bz, that is 


y, =lim,y, u,=lim, 4, etc. (8a) 
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The following notation will also be used 


ol 


PoS Er, Use, w= —*. (8b) 
The slight inconsistency in the notation should be kept in mind: y; is not 
obtained by applying the second limit process to p, but by applying it to p/e 
and afterwards multiplying it by «. If now 7 is again replaced by y/e, ye, etc., 
are functions of &, 7, and ». The flow field which has wy, as stream function is 
the flow field given by boundary-layer theory. [The relation of y; to ue, wz, 
etc., will be discussed below, see (10b) and ff.]. The definition of the boundary- 
layer approximation thus introduced is formally different but actually equi- 
valent to PRANDTL’S definition (cf. also [1], [2], and [3]). It is an approximation 
with respect to a coordinate system ¢. A different coordinate system would in 
general give rise to a different approximation, as discussed in Section 3. 

The functions ye, uz, ve, Me, and p; depend on the two variables € and y/e 
only, and their variation with the latter variable represents the transition 
across the boundary layer. Note that the boundary-layer approximation is 
formally defined wherever the ¢-coordinates are defined, hence a flow field is 
formally given even outside the boundary layer proper. It is this complete 
flow field that will be discussed in the present paper. 

We, Uz, etc. satisfy the Prandtl boundary-layer equations, which actually 
may be obtained by applying the second limit process to equations (1’) or 
(2’). As an example, (2’) become 


oye 9 Ope < 

( ia ’ 0g ~~ 0€ i) =Vele g Ps a a 2 W:, (9a) 
Pak 292 07 We 

an Sere (9b) 


Here the subscript w indicates that the quantity in question has been 


evaluated at the wall: gi“(é) = g°*(¢, 0), etc. 
Note that the relation of y to w and v is unaltered by the second limit 


process since pz, vz, and 7; still satisfy the relations 


ICN: OU ey 10a 

er te Us, kay Ue, ( ) 
and 

Ope tf Ope 10b 

a ae =—4,. (10b) 


As a consequence, (1’) and (2’) remain equivalent after the second limit process 
has been carried out. The flow field given by y; is the same as that given by 
uz and v;. On the other hand, the relation of to this field has been altered. The 
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vorticity of the flow field described by ye is 


~ Wet ae Visi" 3B). 


As is seen by (9b), this is in general different from @;. 

The boundary-layer approximation may also be defined as the solution 
of (9) with certain approximate boundary conditions which will be discussed 
later. For the present, however, it will be regarded as the result of applying the 
second limit process, in the manner described above, to the exact solution of 
the Navier-Stokes equations. 


3. Comparison of Different Boundary-Layer Solutions 


The boundary-layer approximations to an exact solution of the Navier- 
Stokes equations are not uniquely determined, but depend both on the choice 
of dependent variables used to represent the flow field, and on the independent 
variables, i.e., the coordinates. 

Dependent variables. As remarked above, the (w) and the (uw, v) representa- 
tions are equivalent. In general, however, different representations lead to dif- 
ferent approximations. For example, when (&, 7) are orthogonal coordinates, 
the flow field is often represented by means of velocity components (U, V) in 
the sense of vector geometry (cf. [8], p. 101, and especially [7], p. 114). The 
magnitudes of U and V are then equal to the projections of the velocity vector 
on the local coordinate axes. They are connected with the contravariant vector 
densities u, v through the relations 


Zhe Hye. US Ate.) E: (11a) 


where ds? = hj d&? + h? dy? is the square of the line element. Applying the 
second limit process, one obtains 

Up = IhewU,, V=hhy V. ) (11 b) 
where h;,, = h,(&, 0) and V = V/e. Thus the original metrical relations (11 a) 
no longer hold, except at the wall. The flow field corresponding to w; and ¥; will 
then be different from that corresponding to U; and V;. The example of the 
(U, V) representations also shows that, in general, the boundary layer approxi- 
mation does not satisfy the exact continuity equation!). However, the (y) and 
the (u,v) representation have the special property of preserving the exact 
continuity equation under the boundary-layer approximation. In the following, 
only the (y) and the equivalent (u,v) representations will be used. While a 
further investigation of the role of the dependent variables may be of interest, 


1) Cf. Ref. [8], p. 119. 


' 
. 
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4 it will not be carried out in the present paper, which is concerned mainly with 
_ the role of the independent variables. 


i 
‘ 


4 


The flow field will, in general, also depend on the choice of the parameter ¢ 
as function of v. In the following it will therefore be assumed that ¢ = const Vy. 


_ The choice of the constant will not affect the flow field. 


Independent variables. Main correlation theorem. Consider now two boundary- 


4 layer approximations to the same exact solution, one based on the system of 
_ coordinates ¢ = (€, 7) and the other one on another system y = (@, a). Here 
: and o respectively are assumed to be zero at the wall. The two flow fields may 
_ be represented by the stream functions y; and y, respectively. Their relation 
is then described by the following correlation theorem: 


THEOREM 1. If wy; = f(&, 7) is given, then the boundary-layer solution with 
respect to y can be found directly by the substitution formula 


Wy = f(E5 7) (12a) 

where 
E, = lim, € = &(o, 0) , (12b) 
iy = lim, = (52) o (12c) 


The theorem states that the result of applying lim, to y is identical with 
that of first applying lim; and then lim,. 

The evaluations of the limits (12b) and (12c) are easily shown to be correct 
by considering é and 7 as functions of o for a fixed g, and using the fact that 
and o vanish at the wall. Equation (12a) can be obtained directly, if it is 
assumed that y(&, 7, €) is a continuous function of €, 4, and « (in the space of 
those variables), at ¢ = 0. The validity of the theorem can be checked also by 
verifying that p,, as given by (12a), satisfies the appropriate boundary-layer 
equations [cf. (9)] and boundary conditions [cf. (14) and (18) below]. 

The dependence of the boundary-layer solution on the choice of coordinates 
can now be studied with the aid of (12). Since yp is treated as an absolute scalar, 


yz and yw, represent the same flow field if and only if the systems ¢ and y are 
- connected by the relations 


E(o, 0) = &(0, 0), (13a) 
nie) = (Sa)oao see! 

These relations are equivalent to 
E= h(e) , (13') 


n= flo), (13b’) 
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where /, and f, are arbitrary functions. According to (13a’), only the shape of 
the curves € = const, but not their spacing (‘labelling’) has an effect on the 
boundary-layer approximation. Equation (13b’) implies that a relabelling of 
the curves 7 = const, other than 7 > const 7, would alter the flow field 
given by the boundary-layer approximation. On the other hand, certain changes 
of shape of the 7-curves do not change the flow field. 

The flow fields given by y; and yp, will now be compared in the case where ¢ 
and y are arbitrary systems, not necessarily related by (13). 

Conditions at the wall. In general it should be noted that, since y = 4 = 0 
at the wall, the first and second limit processes are identical there. Since the first 
limit is independent of ¢, the same must then be true for the second limit at the 
wall. It is then a direct consequence of the boundary condition (3a) that, at 
the wall, 

u-(, 0) = v,(&, 0) = 0 . (14) 


independently of the system of coordinates. 

Furthermore, we have the important relation that the skin friction obtained 
from the boundary-layer approximation 1s independent of the choice of coordinates. 
This can be seen directly, as follows. Let t denote the exact skin friction (the 
subscript w being omitted from the usual notation for convenience). In 
analogy with (8), define t = t/e, tT; = lim; T, and t; = € T;. Then 


=e lim [e1(E, 8] . (15) 


Since the limit on the right does not depend on the nature of coordinates, 
T, = T, at any point on the wall. 

Whenever the no-slip condition (3a) holds, the exact skin friction is simply 
—j times the vorticity at the wall 


T= Uh Wy > Hoy (5-7), (loa} 


Passing to the limit 
re 99 0? ¢ 
Brie Fi ee peaele. (16) 


Thus, in computing t; from the velocity field given by the boundary-layer 
approximation, one may use the exact relation between stress and rate of 
strain. However, if the no-slip condition does not hold, t; would denote the 
skin friction computed by means of an appropriate approximate relation, 
derived from the exact relation by the passage to the limit. If the exact stress 
‘relation is used, the skin friction may not be invariant when slipping occurs. 
Differences at a distance from the boundary. At the boundary the choice of 
coordinates is irrelevant in the sense discussed above. Within the boundary 
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_ layer proper, that is, within a distance of order ¢ from the wall, the flow associ- 
_ ated with y; and y, should agree up to and excluding terms of order e. This 


L 


_ According to (12), the flow (y,) may be regarded as an image of the flow (wz) 


f 


- follows from Theorem 1. 


However, at a distance from the boundary, », and wy, may differ radically. 
obtained by the following mapping: let a point P go into an image point Q i 
such a way that, if the é- and 7-affices of P are 


&(P)=a, 7(P)=5, (12e) 

_ the point Q is determined by the relations 
| &(Q)=4, &7,(Q) =. (12) 
Then, the stream function has the same numerical values at Q and at P, that is 
Y,(Q) = Ye(P) - (12g) 


Thus, a streamline y, = c will be the direct image of the streamline py; =, 

but, in general, the coordinate curves will not be direct images under this 

mapping. However, if y and ¢ are connected at the wall by the special relations 
0 

Cw = Sw, (52) 4; (12h) 


&, becomes equal to @ and ¢ Ny to a, so that the y-affices of 9 become 
0(Q)=a=€(P), o(Q)=b=7(P). (12}) 


Thus, if (12h) holds, the coordinate curves as well as the streamline curves 
map as direct images. This means that a given streamline, y = c, passes through 
the same coordinate affices in both flows. 

Now, keeping conditions (12h) imposed at the wall (with respect to some 
fixed ¢-system), one may deform the y-coordinates in an otherwise general 
manner. In particular, in any region external to the wall, of the type 7 = > > 0, 
any desired continuous deformation of the y-coordinates can be performed 


without violating (12h). Since the direct-image relations then hold, it is clear 


that, within a region of this type, the family of streamlines may be deformed 


into any desired continuous image of itself by a suitable choice of y. 


Thus, at a distance from the boundary, the flow field obtained from the 
boundary-layer approximation depends almost exclusively on the choice of 
coordinates, rather than on the physical problem. For this reason it is custom- 
ary to disregard the boundary-layer approximation entirely outside the 
boundary layer. The boundary-layer solution is assumed valid as a certain 
approximation up to the edge of the boundary layer, where it is patched on to 
a potential solution. However, the viewpoint taken in this paper is different. 
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One may make positive use of the dependence of y; on the choice of ¢ by search- 
ing for special systems of coordinates, such that y; approaches the exact 
solution y as closely as possible. Then y; would be a valid approximation even 
outside thé boundary-layer region proper. The problem of such ‘optimal’ 
systems of coordinates will be discussed in the following section. 


4. Comparison with Exact Solution. Choice of Optimal Coordinates 


In the preceding section, it was shown how the boundary-layer solution, as 
described for example by w;, varies radically with the choice of ¢. The question 
was raised whether one could turn this fact to good advantage, by choosing € 
in such a way that wp; is in some sense very close to y, even outside the boundary- 
layer region proper. The conditions that such a y; would have to satisfy, need 
to be specified more precisely. It is certainly natural to require that u- and v; 
satisfy the same boundary conditions as w and v at infinity. A stronger require- 
ment would be that, near infinity, py; should agree with y up to a certain 
order in ¢. It is, in fact, possible to require that y and y; agree up to and includ- 
ing first order terms in ¢ everywhere outside the boundary. An equivalent 
requirement is that y; contain, in a sense to be defined, the external flow and 
the flow due to displacement thickness. (The latter flow is actually the first 
order perturbation of y, outside the boundary.) 

Below it will be shown that systems of coordinates € may be found, such 
that yw; satisfies this requirement. Such systems will be called optimal. The 
Discussion (Section 6) will briefly touch upon the question whether other 
requirements on wy; would be preferable. 

Behavior of wz for large. First order approximation to p for n + 0. A funda- 
menta] assumption of boundary-layer theory is that, for a flow property / 

lim (lim, /) = lim (lim, /) (17) 
> 


> w 
provided the limits exist. Applying (17) to « yields 
Uy >Ugy aS 2 >OO. (18) 


Applying (17) to v or v, or to y or p, however, leads in general to trivial results. 
To find the behavior of y; and v; for 7 large, higher order external perturbations 
have to be considered. 


Previously, the external flow had been defined as 
y.=lim,y, 4,=lim,4, etc. (19a) 


We shall now consider this definition to be valid only for 7 +0 and extend 
the limiting functions by continuity to the wall: 


Vow = lim WE 9). 5 ates == LM HAE Te CLC, (19b) 
7 n—>0 


. Vol. V, 1954 The Role of Coordinate Systems in Boundary-Layer Theory 123 


Note that if (19a) is applied even for 7 = 0, y, would have been continuous 
_ there, whereas u, in general would be discontinuous. 

Due to the presence of the boundary layer, it is to be expected that the first 

_ external perturbation is of order e, that is, 


y= y.(P) + € pi(P) + o(e) (20) 


_ where lim, o(¢)/e = 0. Thus the first order external perturbation may be defined 
_ by a repeated application of the first limit process as 


ye = lim, Y= , (7 +0) (21a) 
Vow = eG, 0) = Jim wel, 1) - (21b) 
> : 


Again it is necessary to extend the definition to 7 = 0 by continuity, since 


lim, fu — Yew. =f 


_ whereas Wew, as will be seen by (23), is in general different from zero. 
Now, for 7 small, 


Op n? 
Pe = Usw 4H ( On? Mv 4 SIF Keet 
Hence lim; (y,/é) = ue» 7 and lim; [(y — y,)/e] = pz — ew - Applying (17) 
_ to (yw — »,)/e one then finds that 


Ye— Yew > Pow aS > 0O (22) 
or 
We ~UewN+EWeow for 7 large. (22') 


Since Wr — Uy = Jf (ue — Up) ay it follows from (22) that 


co 


Wore i) (tp — Ue) aN. (23) 


0 


This relation shows how yp, ,, may be evaluated from any given boundary-layer 
approximation. On the other hand, since y, may be defined by (21), the integral 
in (23) must be independent of the choice of ¢. This may also be checked from 
(12). It may be seen that y; satisfies the Laplace equation and hence represents 
potential flow. y, may then be determined if one uses (23) as a boundary con- 
dition, together with the additional condition that y;,/y, > 0 at infinity. It 
follows from the preceding discussion, in particular from (23), that y, repre- 
sents the external perturbation flow due to displacement thickness. y, repre- 
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sents a flow past the body so that the body surface is a streamline (y, = 0) - Wer 
on the other hand, represents a flow that seems to emanate from a source 
distribution at the wall. 

First order approximation to p; for n+ 0. The first limit process may be 
applied to the boundary-layer solution y; in the same manner as was done to y. 
The flow field represented by the limit y;, will be called the external flow 
contained in the boundary-layer solution. Since, for a fixed y, 7 tends to in- 
finity as e tends to zero, one may use equation (22) in evaluating the limit: 


Yee = lim, Pe = Uew 1 - (24a) 


The associated velocity component u;, = 0y;,/0n =u, , may be found 
directly as lim,(#,) for 7 + 0, and extended to the wall by continuity. Thus, 
for a fixed &, the curves of u; versus 7 with e as parameter will have the line 
U = Ue» = const as a limiting curve for ¢ tending to zero. The external flow 
field contained in the boundary-layer approximation has thus a constant -com- 
ponent of velocity. This is of course due to the fact that w- depends on y and € 
only in the combination y/e, so that the limiting value u;, is u-(§, oo). On the 
other hand, the é-component of the exact solution may have a more comphcated 
dependence on 7 and «. Hence, its corresponding limiting curve, u,(€, 7), may 
vary with 7 and need coincide with u-, only at the wall. 

Similarly, one may define the external first perturbation contained in the 
boundary-layer approximation and use (23) in its evaluation: 


Veo = lim, TE Tks = yey. (n +0) (24b) 


Its value for 7 = 0 is defined as the limiting value when 7 +> 0. Since w;, is 
independent of 7 it follows that w:, is zero and vz, = — Oy, ,,/0f = v,. Hence 
ve, and v;, regarded as functions of 7, must coincide for 7 = 0, but as 7 varies, 
the former must remain constant, whereas the latter may vary. 

Optimal coordinates. The boundary-layer solution will be said to contain 
the external flow of the exact solution, whenever 


Wee = Ye i@n, Usw 1] = We (25a) 
and to contain the flow due to displacement thickness, whenever 
Pre = We Le, Yow = Yo- (25b) 


A coordinate system will be called optimal if both these equations are satis- 
fied. Assuming that y, and y; are given, let y = (@,¢) be such an optimal 
system. It can be seen immediately that (25b) is satisfied if and only if y} is 
a function of @ only. In other words, it is sufficient (and necessary) to choose 


; Vol. V, 1954 The Role of Coordinate Systems in Boundary-Layer Theory 125 


_ the curves @ = const to coincide with streamlines of the flow due to displacement 
_ thickness. (25a) is satisfied if and only if u, = dy,/0n = 1,4, so that u, is a 
- function of g only. This can be accomplished by taking o = we. Most generally, 
one may choose o = f(g) y, where f(g) is an arbitrary function of 9. ben 
Ue, = 1/f(0) = Ue. Hence, the following theorem has been proved. 


_ THEOREM 2. The coordinate system ¢ = (&, n) is a particular optimal system, 
in the sense that (25a) and (25b) are satisfied, if 


E=w,, 1, : (26a) 


re other system, ¥ = (0, 0), 1s optimal if and only if it is related to the above 
system by 
e=A), o=f(6) 4 (26 b) 


_ where }, and f, are arbitrary functions. The flow field given by the boundary-layer 
approximation ts the same for all optimal systems, but will be different if any 
_ other system 1s chosen. 


The last statement follows directly from a comparison of (26b) and (13’). 
It will now be shown that the optimal boundary-layer solution, as formally 
_ defined by the requirement (25), also has the other properties discussed at the 
beginning of this section. First it will be checked that it satisfies all the bound- 
ary conditions imposed on the exact solution. The conditions at the wall (3a) 
are satisfied according to (14). At infinity y approaches a potential solution 
which actually is y,, whereas y; approaches w;,. If now y, = pz, it follows that 
yr satisfies the correct boundary condition. Note that, according to the same 
boundary condition, w; should approach u,(&, co). On the other hand, for any 
boundary-layer solution, u;(£, 0c) = u,(&, 0). To satisfy the correct boundary 
condition at infinity for ug one should then make u,(&, 0) = u,(&, co). An easy 
way of achieving this is to choose 7 =, as in (26a), since in this case 
_u, = const = 1. 

So far only the relation y;, = y, has been used. If in addition y;, = y, 
-then y and y; agree to first order at infinity and actually in any region which 
excludes the boundary. This follows directly from the fact that for 7 + 0, 


y=y.tey' + o(e) (20) 
and also 
Yr = Vro + E We + O(€) - 


A further evaluation of the result will be attempted in Section 6. 
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Construction of optimal approximation. In general, the optimal system cannot 
be determined a priori. However, one may proceed.as follows. First the external 
flow y, is found by solving the proper potential flow problem. Then a boundary- 
layer solution w; with respect to any convenient system of coordinates is found 
by solving the boundary-layer equations (9), or equivalent equations, using 
(14) and (18) as boundary conditions. By evaluating the integral in (23) one 
finds y;,,, and then the flow due to displacement thickness by solving LAPLACE’s 
equation for y;. y, and y; being known, an optimal system is directly given by 
(26). Finally, to find w, it is not necessary to solve the corresponding equation. 
y, may be obtained from y; by a direct substitution in accordance with (12). 


5. Examples 
Example 1. Flow Towards a Flat Plate 


The stream function 


co p=V9U x1Q), (27a) 
where y = VU|y y, and 
(Meh fen (gah (27b) 
(0) = #0) =0; fee)=1, feo) ~¥—e, (27¢) 
e~ (1-1 Olay, (274) 


represents an exact solution of the Navier-Stokes equations for y = 0, satis- 
fying the conditions y = dy/dy = 0 at y=0, y~Uxy—cxVyU as yoo 
(cf. previous comments and equation (3c)]. Thus y represents a flow normal to 
an infinite wall, as in Figure 1 (see [1] and [2]). 


Y 


Figure 1 


Holding U, x, and y fixed, and letting » +0, we obtain the external flow. 


Thus, : 
y. = lim, ou x2 y)| =Uxy.. (28) 
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aoe ee 


: Choosing the boundary-layer thickness parameter as ¢ =// r/U, and regard- 
ing U as constant under the limit process, the boundary-layer approximation 
with respect to the (x, y)-coordinates is found by applying the lim, process. 
Thus, 


ey ar 


Pm TAO 


y, = lim, © = lim, U x f(y) = U x f@) (29a) 
or 
q ema Vy ee 
a Aes ee vUx fly). (29b) 
; Finally, the flow due to displacement thickness is found as: 
q 2 lin (eee ok ”) arts ee (30) 


_ Thus, the optimal &-coordinates are € = /,(w,), or, since the choice of f, is in- 
_ consequential, € = x. The optimal 7-coordinates are 7 = f,(&) y,. Since the 
_ choice of f, does not affect the result, we can take /,(€) = 1/(U &), that is, 
_»= yy. Thus, in the present case, the (x, y)-coordinates lead to the optimal 
_ boundary-layer approximation, which, as may be seen by comparing y, with y, 
happens to satisfy the full Navier-Stokes equation. 


e Example 2. Oseen Flow Past a Semi-Infimte Flat Plate 


hae. ee 


In this example the Oseen equations will be regarded as a model for the 
_ Navier-Stokes equations, rather than as an approximation. In other words, 
the Oseen equations will be taken as the fundamental ‘exact’ equations and 
_we shall be concerned with approximations to these equations. The results 
_ of the present paper apply, of course, equally well to the Oseen equations as 
_ to the Navier-Stokes equations. 

The Oseen equation for y is: 


(xy? U =) ty =0. (31a) 


- The boundary conditions for a flat plate are 


0 
Vary, = 0 ade Oe ge > 10) | (31b) 
: as Oy x 9 As bey eae OO, OT as’ x —0o . (31c) 
f Oy Ox 


The stream function y, is obtained from the equation resulting by the first 
limit process, which is 
; 0 
OF y>y.=0, (32a) 
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and from the boundary conditions: 


0 Oy 
ee Ol Sat she Oa ey, Sapo Us SF aakas as | vy| +00, x > —oo 
(32b) 


(cf. Section 2), corresponding to the assumption of no separation. Furthermore, 
in view of symmetry of the viscous problem about the x-axis, it is expected 
that the circulation is zero. Thus, 


ye =Uy. (32c) 
Next, choosing ¢ = V»/U, (31a) is written in the form: 
Oe 0? 0 0? gq 02 Nee 
iCex. = os poe oat) ¥ se e233 


Holding x, y (and nominally, U) constant, and letting «>0, yp > y,, we 
obtain the following boundary-layer equations: 


02 0\ Op, — ie 
ser rd ae i33B) 
and the conditions: 

p= GE =0 at p=0,x>0. '(33c) 


Furthermore, at infinity, the boundary conditions are: 
~~ UY + Yow (33d) 


where yw, 1s unknown. Hence, the boundary-layer approximation with respect 
to rectangular coordinates is: 


nee =| 9 7 1 (* | 
= 2U\V x ae erf Ae a Ten ae i) (33e) 


From equation (33e) we now find 


Pew = — : (34a) 


Putting y= y, + ey’ in equation (31a) and letting « > 0, the perturbation 
equation is 


0 ' 
- V2w,=0, (34b) 
or, since wy represents an irrotational flow as « > --oo , 


Vey, =0. (34) 


eee ee ee 


Bat a AN i Dl a nl | ds | 


PEt IN NRE ee 


Dente ce), ee 
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The solution corresponding to a distribution of fundamental sources at the wall, 
which satisfies (34c) and (34a) is: 


, + *¥ 
v= 2 |/ 7. (34d) 


Since y, is constant along confocal parabolas, the optimal &-coordinates 
may conveniently be taken as: 


&= REG 5 C=Vat+iy, 0 sSarg(4 +41 y)< 22. (35a) 


The optimal 4-coordinates are of the form 7 = f(€) y, or, choosing /(&) = 1/(2 &), 
Wwe may take 


y= IC. (35b) 


The classical approximation is now transformed to optimal coordinates 
using the transformation (12). The quantities x;, yz, are, from (35): 


BF), (36a) 


S 


Ye = 287. (36b) 


_ Hence, the optimal approximation is: 


y= 208th et] + T= lasts | (36) 


or 
1 


Va 


While it can be seen that the classical approximation has a non-vanishing 


ye = 2V 0 v8 [ert + (er —1)}. (36d) 


_ v-component of velocity at infinity and does not satisfy all the boundary con- 
_ ditions for y, and is further characterized by all the ‘anomalies’ of the Blasius 


solution, it can be verified that the optimal boundary-layer approximation 


_ satisfies exactly the Oseen equation and all the boundary conditions (31)’). 


Example 3. Flow A gainst a Wedge 


For this class of flows, the exact solution of the Navier-Stokes equation is 


_ known only when the wedge angle is 180° (cf. Example 1). y, is therefore found 
_ by solving the corresponding potential problem, which gives 


2—-—% | -jnal/(2-0 2/(2—a 
yes Tm [US einen) (eaypt-2)), (37) 


1) The solution of OsrEN’s equations is due to Lewis and CARRIER [9] and the phenomenon 
that the boundary-layer approximation in parabolic coordinates is exact, has first been noticed 
by VAN Dyke [4]. 
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where z* = x* + 7 y* is the complex variable in the physical plane, U is the 
speed of the external flow at z* = —1, and wa is the apex angle (cf. Figure 2). 


* 


F 


Figure 2 


The boundary-layer solution for the wedge has been investigated, among 
other authors, by H. WEyL [1], who uses ‘streamline-potential’ coordinates 
as a basis for the approximation, that is, the coordinates 


x= ag es Oe , (38) 


where q, is the harmonic conjugate of y,, normalized to be zero at the apex. We 
then have y, = U y, or, u, = U = u,y. The boundary-layer equations based on 
the streamline-potential coordinates admit a solution of the ‘similarity’ type, 

namely, of the form | 


@.= UV 2fQ). A= A (39a) 
where y = yVU/p and 


1,(0) = 7,(0) =0, > f,(eco)=1. (39b) 


Furthermore, as 1 > 0, f,(A) ~A —c, where cy = f{t — fi,(A)] da. 
0 


THUS: oy = =U Cy Vx, and, therefore, solving LAPLACE’S equation, we 
obtain : 


w,= —Uc,&, where E=ReVx+tiy. (40) 


It is thus seen that the streamline-potential coordinates (38) are not optimal. 
They have the limited advantage that y, is contained in the boundary-layer 
approximation, but y, is not contained, and all the singularities of the Blasius 
approximation are present. 

The optimal €-coordinates may be taken as: 


E=Re€ (41a) 
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and the optimal 7-coordinates may be taken as: 


RAYE Oo NET NS. 


ie Ina (41b) 


where € = Vz. This gives a conformal set of coordinates, which may be called 
_ streamline-parabolic. The optimal coordinates in Examples 1 and 2 are special 
_ cases of this system. 

) We have x; = €?, yp = 2&7, and so the optimal boundary-layer approxi- 
_ mation is: 

¢ Y= UET,(2 7%) , (42) 


_ where /, is the same function as that in equation (39). 

: Now, as 7 > 00, pp ~2UEn —eUc,é=y, + € yy, and further, the solu- 
tion extends over the entire plane and is free from anomalies. Unfortunately, 
however, the Navier-Stokes equation is not satisfied by the optimal approxi- 
mation, except if « = 1, corresponding to Example 1. 

A special case is x = 0, which represents flow past a semi-infinite flat plate, 
in other words, the non-linear version of Example 2. In this case parabolic 

_ coordinates are optimal. They have been used for this case in [10], [11], and [12] 
where, however, the inclusion of flow due to displacement thickness is not 
commented upon. 

It should be pointed out that, in general, the optimal coordinates cannot be 

{ taken as conformal. This is evident, for example, if the streamlines y, = const 
meet the wall at angles other than normal. An example of this is flow past a 
finite flat plate. 


6. Discussion 


The significance of Theorem 2 will now be discussed. 

First, it is obvious from (15) that if one is interested only in the first 
approximation to the skin friction, then the choice of coordinates is irrelevant. 
The purpose of choosing an optimal system is rather to obtain a good picture 
of the complete flow field. The problem studied here, is, vaguely speaking, that 

of obtaining as much information as possible within the limits of boundary- 
layer theory. It is reasonable to require this information, before proceeding to 

| approximations which are definitely beyond boundary-layer theory, for 
example, in considering higher order approximations to the skin friction. 

The flow field of the optimal solution will therefore be compared with a 
flow field obtained by standard methods, that is with the composite flow 
field defined as follows: A certain definite edge, terminating the boundary- 
layer region, is introduced. Within the boundary-layer region, the flow field is 
taken from a non-optimal boundary-layer solution ; outside the edge, a poten- 
tial flow, consisting of the external flow and flow due to displacement thickness, 
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is assumed. The viscous and non-viscous solutions are patched together along 
the edge. 

First a general observation can be made. The unpatched solution, that is, 
the optimal boundary-layer solution, is conceptually clearer. It is free from ~ 
artificial singularities and therefore easier to handle analytically. The question 
of where and how to patch and the problems associated with this are avoided. 

In order to make a more detailed comparison it is convenient to consider 
separately the potential part of the flow and the viscous shear waves. These 
terms are easily defined in the case of Oseen flow, since in that case the velo- 
city field separates into a potential and a transversal wave in such a way that 
all pressure forces are associated with the former .and the viscous shear only 
with the latter. A suitable definition for the non-linear case will be discussed 
later. For the time being we may consider the vorticity as representing the 
viscous shear waves. Thus consider first the vorticity. In the patched solution 
it is assumed to be zero outside the boundary layer. In the optimal solution it 
vanishes continuously away from the boundary, but very rapidly (exponen- 
tially in Example 2), so that the actual numerical correction is small. However, 
in the neighborhood of the leading edge, the difference may be large. Consider, 
for example, the flat plate (Examples 2 and 3). If rectangular coordinates are 
used, the entire line x = 0 is singular. The vorticity tends to infinity as x de- 
creases to zero for any value of y. Furthermore, no vorticity is obtained in the 
region upstream of the leading edge. These defects are remedied by the use 
of parabolic coordinates. The artificial singularity along the line x = 0 is re- 
moved, and a certain upstream vorticity is obtained. In the Oseen case, the 
optimal solution is, in fact, exact and thus holds in particular near the origin. 
In the non-linear case, the optimal solution is not correct near the origin, but, 
qualitatively, comes much closer to the anticipated actual behavior of the exact 
solution then the solution in rectangular coordinates. It follows from the 
elliptic nature of the Navier-Stokes equation that an upstream spreading of 
vorticity actually must occur. Similar remarks apply to the wedge (Example 3), 
in comparing the solutions in streamline-potential and streamline-parabolic 
coordinates. Inside the boundary-layer region no improvement in the values 
of the vorticity may be guaranteed by the use of optimal coordinates. This will 
be commented on later. 

Next consider the potential part of the flow. Outside the boundary layer 
proper the patched flow field is, by construction, simply the potential part of 
the optimal solution. In considering the flow inside the boundary-layer region 
it should be remembered that any boundary layer solution contains part of the 
external flow and flow due to displacement thickness. In the terminology of 
Section 4, the relevant comparison is thus between yp, + € y, and Pte + E Vre 
the former being the correct value (to the order e). In particular it is of interest 
to compare the corresponding first derivatives since these represent the 


4 
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velocity fields. y;, and its first derivatives have the correct values at the wall, 
but their normal derivatives of higher order are different. Hence, within the 
boundary-layer region (of thickness Vv) ue, may show a mistake of order Vy. 
Similarly, only the values of y:, are correct at the wall, while the normal deriv- 
atives are not generally correct. Thus, another error of order Vy generally 
appears in w;,. These errors cause difficulty in connection with the patching 
et the edge of the boundary layer. 

It should be remembered, however, that the criterion of optimality adopted 


here is still somewhat arbitrary. A more systematic approach would be to use 


boundary-layer theory as a basis for finding the beginning of an expansion of 


_-the form 


y~ polP, ») + eywl(P, ») ++: (43) 


where each y, is of order unity uniformly in the complete flow field. (Certain 
singular points, like the leading edge, would require special consideration.) To 
the author’s knowledge, the first systematic discussion of an asymptotic ex- 
pansion of the type (43) for the case of singular perturbation problems was 
given by G. E. Latra in [3]. Considering the optimal solution in this light, we 
see that it is valid uniformly up to and including terms of order ¢ in any region 
which does not include the boundary, since here only the external terms 


We + € y, are involved. On the other hand, if the region includes the boundary, 


the optimal solution is no¢ necessarily accurate enough to include terms of 
order e since terms of order ¢ are in general omitted fromthe boundary-layer 
equations. Actually, the improvement of the boundary-layer solution discussed 
in this paper was based mainly on a consideration of the potential part of the 
flow field, even if some incidential improvement of the viscous shear waves 
resulted. To investigate the matter further it then seems natural to treat the 
viscous shear waves and the potential part of the flow field separately. 

These terms, used previously in a somewhat vague fashion, may be defined 
as follows. The potential part is the part obtained by a repeated application of 
the first limit process, namely 


" . 1 7 
y,+ey,ter?y, +-::, where y, = lim, |< (py — y.—ey,)|, etc. 


The viscous shear waves would be the remainder, namely the difference between 
y and its potential part!). Approximations to the potential part may be found 


' by successively solving LAPLACE’s equation with appropriate boundary con- 


1) The series representing the potential part may be divergent. In this case, the procedure 
would only amount to a splitting of the asymptotic expansion (43), term by term, into potential 
and viscous parts. On the other hand, it may be possible to define the potential part exactly, by 
considering the behavior at » fixed and 7 (= x? — y2) large, that is, by subtracting from yp 
harmonic terms of all finite orders in (74) as y > 00. 
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ditions (which for the higher terms are related to the shear waves). Boundary- 
layer theory may then be applied for the purpose of finding an approximation 
to the viscous shear waves. This procedure, in a sense, represents patching 
at the wall, rather than at some edge of the boundary layer. An equivalent 
procedure is actually used in [3]. If this method is followed the choice of coor- 
dinates should be guided by a study of the shear waves alone. In general, 
different coordinates may then appear optimal. It should be noted, however, 
that for Example 2 and certain special cases of Example 3, Latta was led 
to the introduction of parabolic and streamline-parabolic coordinates from a 
consideration of a series expansion of the form (43). Thus in these special cases 
different methods yield the same results. A closer integration of the methods of 
[3] and the methods used in this paper would undoubtedly prove very fruitful. 

Finally, a brief comment will be added about the case when separation 
occurs at the body. The use of boundary-layer theory in this case is hampered 
by the fundamental difficulty that the relevant external flow y, is different 
from the potential unseparated flow past the body. Then it is actually not 
correct to determine the point of separation from the pressure distribution 
corresponding to the latter flow. In principle, however, if y, is known or 
assumed, boundary-layer solution may still be regarded as a singular pertur- 
bation of y, for y > 0. Asemi-optimal solution, containing at least y,, may then 
be found by choosing suitable 7-coordinates [cf. (26)]. 

A similar problem arises in connection with y; when the flow separates at 
the body. The external disturbance is now influenced by the behavior of the 
separated region. It is not even clear, a priori, if the first external perturbation 
would be of order ¢. However, if the disturbance is of this order, it may still be 
included, in principle, in the optimal boundary-layer solution by a suitable 
choice of &-coordinates. 

Furthermore, it is believed that a complicated singularity may occur at the 
point of separation, at which the v-velocity becomes infinite [13]. It may be 
possible that this infinity can be removed by an adequate choice of the é-coor- 
dinates, in analogy with the removal of the infinity along x = 0 in Example 2. 
This possibility, however, has not been investigated. 


REFERENCES 


[1] H. WeEyL, On the Differential Equations of the Simplest Boundary Layer 
Problems, Ann. Math. 43, No. 2, 381-407 (April 1942), 

[2] R. von MisEs and K. O. FRIEDRICHS, Fluid Dynamics, Chapter IV (Brown 
University | [Notes], Summer, 1941). 

[3] a 1D aie Singular Perturbation Problems, Thesis, Calif. Inst. of Technol. 

O50). 
[4] M. D. VAN Dyke, Private communication. 
[5] P. A. LaGrRsTRom, Private communication. 


i a 


, 
r 
2 


iN) 


TT} 
ae 


“a 


Vol. V, 1954 The Role of Coordinate Systems in Boundary-Layer Theory 135 


[6] P. A. LaGERstRom and S. Kapiun, The Role of the Coordinate System in 
Boundary Layer Theory. VIII International Congress for Theoretical and 
Applied Mechanics (Istanbul, 1952). 

[7] L. Brittouin, Les Tenseurs (Dover Publ., New York, 1946). 

[3}"S: ees ed., Modern Developments in Fluid Dynamics, Vol. I (Oxford, 
1938). 

[9] J. A. Lewis and G. F. Carrier, Some Remarks on the Flat Plate Boundary 
Layer, Quart. appl. Math. 7 (1949). 

[10} Akademia Nauk S.S.S.R., N. E. Kocurn, Sobvanie Sochinenii, Vol. II, 
Collected Works (published by U.S.S.R. Academy of Sciences, Moscow and 
Leningrad, 1949). 


[11] H. L. ALDEN, Second Approximation to the Laminar Boundary Layer Flow 


Over a Flat Plate, J. Math. Phys. 27, No. 2 (1948). 
[12] G. F. Carrier, and C. C. Lin, On the Nature of the Boundary Layer Near 
the Leading Edge of a Flat Plate, Quart. appl. Math. 6 (1948). 


_ {13] S. GoLpsTEIn, On Laminar Boundary Layer Flow Near a Position of Sepa- 


> 


vation, Quart. J. Mech. appl. Math. 7, Part I (March 1948). 


Zusammenfassung 


Die Grenzschichtnaherung zu einem gegebenen Str6émungsproblem ist nicht 


invariant, das heisst, sie hangt von dem Koordinatensystem ab, in dem der 
_ Naherungsprozess durchgefiihrt wird. Zwar ist die von der Fliissigkeit ausgetibte 


Reibungskraft von der Wahl des Koordinatensystems unabhangig, aber das 
Strémungsfeld selbst hangt im allgemeinen vom Koordinatensystem ab. Theo- 


rem 1 besagt, dass die Grenzschichtlésung eines Koordinatensystems durch eine 
- einfache Substitution in die Lésung irgendeines anderen Systems tbergeftihrt 


werden kann. Nur wenn die zwei Koordinatensysteme durch eine ganz spezielle 


_ Beziehung verbunden sind, bleibt das entsprechende Strémungsfeld dabei un- 


wm 


verandert. Diese Abhangigkeit des Str6mungsfeldes von der Wahl des Koordi- 
natensystems wird benutzt, um Theorem 2 zu beweisen, das besagt, dass Koor- 
dinatensysteme gefunden werden kénnen, die in dem folgenden Sinne optimal 
sind: Die entsprechende Grenzschichtlésung ist im ganzen Raum, auch ausser- 
halb der eigentlichen Grenzschicht, als Naherung giiltig. Sie enthalt die aussere 
Potentialstromung und auch die Strémung, die von der Verdrangungswirkung 
der Grenzschicht herriihrt. Durch dieses Theorem wird das Problem hinfallig, wie 
man die aussere Potentiallésung und die innere Grenzschichtlésung an irgend- 
einem Rande der Grenzschicht zusammenfiigt. Die vorliegende Arbeit beschrankt 


sich auf Str6mungen ohne Ablésung, weil die Grenzschichttheorie im Ablésungs- 


falle mit prinzipiellen Schwierigkeiten behaftet ist. Es wird ausserdem voraus- 


. gesetzt, dass die Strémung zweidimensional, stationar und inkompressibel ist. 


Die letzteren Beschrankungen sind weniger wesentlich und werden in einer dem- 
nachst erscheinenden Arbeit diskutiert. 


(Received: June 7, 1953.) 


136 7s ZAMP 


Contribution a l’étude des compresseurs axiaux supersoniques 


Par PIERRE SCHWAAR, Zurich?) 


I. Introduction 
Définitions et généralités 


Dans un compresseur axial ordinaire, les transformations énergétiques su- 
bies par le fluide véhiculé dans l’ailetage s’effectuent en régime d’écoulement 
subsonique. En particulier, l’énergie de pression est engendrée par ralentisse-_ 
ment continu de l’écoulement relatif dans les aubages. Au sens actuel du terme, 
un compresseur axial est dit «supersonique» lorsque son ailetage est le siége 
d’ondes de choc au moyen desquelles l’énergie cinétique communiquée au fluide 
par les aubages moteurs est partiellement ou entiérement transformée en énergie 
de pression. Considérons un étage de compresseur axial supersonique: les ondes~ 
de choc, frontales ou obliques, peuvent se produire dans l’aubage moteur, dans 
l’aubage directeur ou dans les deux simultanément. Quand elles sont frontales, 
écoulement relatif aval est subsonique; quand elles sont obliques, cet écoule- 
ment peut étre subsonique ou supersonique. Un aubage est dit «transsonique » 
lorsque l’écoulement relatif du fluide y passe du régime supersonique au régime 
subsonique par l’entremise d’une ou de plusieurs ondes de choc”). Un aubage 
ou cette transition s’effectue par ralentissement continu de l’écoulement ne 
présente actuellement aucun intérét pratique. Enfin, un aubage est dit «sub- 
sonique», respectivement «supersonique» lorsque |’écoulement relatif du fluide 
y est enti¢rement subsonique, respectivement supersonique®), 

On distingue ainsi toute une gamme de compresseurs axiaux supersoniques, 
dont trois du type a ondes de choc frontales semblent retenir plus spécialement 
V’attention des turbomécaniciens: 


a) Aubage directeur subsonique, aubage moteur transsonique; 
6) aubage directeur transsonique, aubage moteur subsonique ; 
c) aubage directeur transsonique, aubage moteur supersonique. 


L’étude des compresseurs axiaux supersoniques se justifie par le rendement 
élevé de la compression par ondes de choc (figure 1). Mais la réalisation d’une 
telle machine pose quelques problémes ardus. Celui que nous traiterons dans 


1) Société Ardag S.A., Développement et Recherches appliquées. 
*) Cette définition ne s’applique qu’a un aubage de compresseur. L’aubage directeur d’une 
turbine de LAvAL est un aubage transsonique. 


3) Les définitions de cet alinéa n’ont évidemment rien d’absolu. Elles correspondent a l'état 
actuel de la technique. 
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cette note est le premier a résoudre: c’est le probléme fondamental de l’écoule- 
‘ment du fluide dans l’ailetage. Nous limiterons notre étude au cas de la com- 
“pression par ondes de choc stationnaires frontales. 

Le probléme aéro-thermodynamique que nous allons aborder différe de celui 
qu’on rencontre dans le dessin d’aubage du compresseur axial subsonique en 
ce que la présence d’ondes de choc dans l’ailetage introduit un mécanisme 
{ 
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Figure 1 
Rendement de la compression par onde de choc (k = 1,4). 


particulier qui impose a l’écoulement de révolution des conditions cinématiques 
‘et dynamiques bien déterminées et trés restrictives. 

Il s’agit tout d’abord de réaliser l’onde de choc frontale. Nous supposerons 
a cet effet que la veine fluide est limitée par des parois cylindriques coaxiales 
et que l’écoulement giratoire qui aborde l’onde de choc se trouve en équilibre 
radial, c’est-a-dire qu’il s’effectue sur des cylindres coaxiaux a l'amont immeé- 
diat de l’onde de choc, ce qui n’introduit pratiquement aucune restriction im- 
portante. On voit alors que chaque surface de choc frontale, limitée par 
V’intrados et l’extrados de deux ailettes adjacentes et par les parois cylindriques 
intérieure et extérieure, est une portion d’hélicoide radiale 4 pas constant. En 
effet, considérons un rayon traversant l’écoulement a l’endroit du choc. Dans 
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V’hypothése que nous venons d’admettre, chaque vecteur vitesse qui le coupe 
lui est perpendiculaire. La surface de choc frontale devra donc contenir ce 
rayon. Mais l’écoulement giratoire présentant une symétrie de révolution’), la 
surface de choc est une hélicoide radiale 4 pas constant. Réciproquement, il est 
évident que c’est la seule surface dont les normales définissent un champ de 
directions qui coincident avec celles des vecteurs vitesses d’un écoulement de 
révolution en équilibre radial, écoulement dont la fonction de gauchissement 
est ainsi déterminée d’avance. En fait, l’équation de l’hélicoide radiale en 
coordonnées cylindriques étant: 


h 
Z=mp avec mM=—>—, 


ot h est le pas, la pente principale de l’hélicoide est donnée par 


Roe 
84 dp 
La fonction de gauchissement de l’écoulement giratoire est donc définie au 
signe prés par 


tga(8) =o =a ter. (1) 


On retrouve ainsi la distribution radiale des angles des écoulements gira- 
toires classiques du type a potentiel (c,7—= const, c,»,— const). C’est une 
restriction importante. 

Considérons alors les conditions dynamiques qui régissent 1’évolution trans- 
sonique envisagée: tout écoulement vérifiant la relation (1) et la condition 
d’équilibre radial en amont de l’onde de choc (en particulier l’écoulement 
irrotationnel), ne sera pas en équilibre radial 4 son aval. Il se produira done 
généralement un brusque déséquilibre a la traversée de l’onde de choc, avec 
décollement éventuel de la veine et toutes les conséquences que cela comporte. 
Dans le cas de l’aubage moteur transsonique, c’est non seulement la transfor- 
mation é€nergétique par choc qui s’en trouvera affectée, mais l’énergie elle- 
méme communiquée au fluide par l’aubage. 

La réalisation d’un compresseur axial supersonique 4 ondes de choc fron- 
tales et rendement élevé dépend ainsi des conditions suivantes: 

Existe-t-il un écoulement giratoire supersonique qui vérifie la relation (1) 
et la condition d’équilibre radial des deux cétés de l’onde de choc ? Cas échéant, 
le nombre de Macu local de cet écoulement variera suivant une loi bien déter- 
minée le long du rayon avant (et aprés) le choc. A quelle distribution radiale 
d’énergie d’arrét cette loi correspond-elle? Peut-on réaliser cette distribution 
d’énergie par l’entremise d’un ou de plusieurs aubages moteurs (nécessairement 


1) C’est une hypothése fréquente et difficile A éviter. 
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‘subsoniques ou supersoniques au sens défini plus haut) qui effectueraient sur 
le fluide un travail adéquat fonction du rayon? Enfin, la condition d’équilibre 
‘radial n’impliquant pas la stabilité radiale, les divers écoulements envisagés 
Templissent-ils cette derniére condition? Remarquons encore qu’au cas ot le 
choc se produit dans l’aubage moteur, c’est au moins un étage préparatoire 
entier qui précédera les étages transsoniques. 

_ Lacompression transsonique par ondes de choc frontales pose encore, pour 
tous les types d’aubages utilisables 4 cette fin, un probléme particulier sur le- 
quel nous reviendrons dans la conclusion de cette note: le probléme essentiel 
de la stabilisation de l’onde de choc dans I|’aubage. 

Passant a l’étude analytique des questions que nous venons de formuler, 
on pourrait se proposer de traiter simultanément les cas de l’aubage moteur 
et directeur transsoniques. On passerait du premier au second en annulant la 
vitesse de rotation et en remplacant les symboles des grandeurs relatives par 
ceux des grandeurs absolues qui caractérisent l’écoulement. Malheureusement, 
la géométrie des triangles de vitesses fait intervenir des relations peu maniables 
dans l’expression de la condition d’équilibre radial, de sorte que le cas de 
lYaubage moteur transsonique ne donne pas lieu a une solution explicite. 

Nous étudierons donc premiérement le cas ot l’onde de choc se produit 
dans l’aubage directeur. Sans nous prononcer ici sur les avantages et incon- 
vénients probables des différents types de compresseurs [1], [2]1), disons seule- 
ment que c’est le type c) 4 aubage moteur supersonique que nous avons plus 
particuliérement en vue. Nous reviendrons plus tard sur le cas de l’aubage 
moteur transsonique. 


Les équations et le rendement du choc frontal 


Nous rappellerons tout d’abord briévement les équations qui rendent compte 
du phénoméne du choc frontal. Ce sont: l’équation de continuité: 


é%),  Q 


S»5 


ag 
Véquation des quantités de mouvement: 
precr=p+oct (3) 
et l’équation de l’énergie: | 


r ae 
Sed it (oe inthaden 5 (4) 


1) Les chiffres entre crochets renvoient 4 la Bibliographie, page 150. g3 
2) La signification des symboles est claire ou est indiquée dans le texte. Les lettres qui désignent 
les grandeurs d’état thermodynamiques derriére l’onde de choc sont pourvues d’un accent cir- 


conflexe. 
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En éliminant alternativement deux des six variables #, 0, c, B, 0, C, on obtient 
une série d’équations du deuxiéme degré a deux groupes de solutions: les unes 
sont triviales: 

p=, 0=6, c= (iln’ya pas de choc); 
les autres constituent les équations du choc frontal. Nous les écrirons sous forme 
sans dimension, en introduisant a cet effet les nombres de MACH locaux: 


out k=c,/c, est l’exposant adiabatique, que nous admettrons évidemment 
constant. On obtient ainsi [3], [4]: 


glides it pO Gy Ca 6 
in [Leo icles 1], (HuGonIoT) (6) 

tai Pie 2k ie 
Sieg hast a argo ane (7) 

@ 2 kit 
eo eae? ick Ge (8) 

En outre: 
T° Sa + i ae he? le (0) 
[Sid (k + 1)? M2 : 

: ae ay es 

May Ie a (10) 


On définit enfin le rendement adiabatique de la compression par onde de choc: 


travail de oes oe) paces de P ap 


UCP ach oe ry aR neern ORE, Am 


Les relations de DE SAINT-VENANT et d’HuconiotT permettent d’écrire [5]: 
pprtaiia 1] — 1 P 1| 
1 See Se, : (11) 


La figure 1 donne a titre de rappel le graphique de 7, en fonction du taux de 
compression statique P = £/p pour k = 1,4. 


II. L’écoulement transsonique dans l’aubage directeur 
La condition d’équilibre radial s’exprime par: 


dp Ca 
ay ee oa (12) 
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a 


_ Nous transformerons cette relation en introduisant l’angle « de l’écoulement 


nied) \ \ dees a 


avec la direction périphérique et les variables sans dimension M et y = 7/r;,, 
ou 7; désigne un rayon d’aubage de référence que nous prendrons égal au rayon 
intérieur. Un calcul élémentaire donne: 


d cos? 
ce eR (12a) 


Logarithmons et dérivons (7) par rapport a y et introduisons la condition 


d’équilibre radial (12a) devant et derriére l’onde de choc!): 


ere 2 2k dM? 
bM?2 COS*H g COS" a : 
ee ee TSC ap 


Remplacons dans cette relation M2 par sa valeur (10). On obtient, toutes 
_ réductions faites, |’équation différentielle : 


dy dM? 


Portons enfin dans (13) la condition (1): 
2 i 
tga=Cy, c’est-a-dire cos?a = Te Cry: (1a) 
Aprés décomposition en éléments simples, l’équation différentielle prend la 
forme: 


7 eC oe k+1 


dy C2 y dy 1 dM? dM? 2 
e Fe =e Le oF o at) 


dont l’intégration est élémentaire: 
y Mbt BOS ah 


GT ear oyaacy |e laa pae 1 | 05) 


c’est-a-dire 


Reo 
Tee 
M*(y) =1+ M? + ihe | a y VVRe1 4 (16) 
Teh ezor eae 


C’est la fonction M2(y) cherchée qui détermine avec la relation (1a) l’écoule- 
ment giratoire supersonique qui assure une onde de choc frontale et satisfait 
aux conditions d’équilibre radial prescrites. Pour y oo, M? tend vers une 
limite bien définie M2_,,, =1 fonction de M? = 1 et de C. L’écoulement est 


you = 


1) Le probléme de l’équilibre radial derriére l’onde de choc a été mentionné par W.A. Logs [1]. 
A notre connaissance, aucune étude analytique des influences radiales n’a été publiée jusqu’a 
ce jour. 
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donc supersonique dans tout l’espace annulaire 1 S y Soo. Fait trés impor- 
tant, car il montre, contrairement a certaines prévisions [1], qu’il n’existe en 
principe aucune limite aéro-thermodynamique a la valeur du rapport carac- 
téristique v = 7;,/r7, d’un aubage directeur transsonique. La hauteur relative 
des ailettes et, par conséquent, le débit spécifique!) de l’étage seront cependant 
limités par les contraintes mécaniques de l’aubage moteur?). 

Nous comparerons la distribution M?(y) donnée par (16) a celle d’un écoule- 
ment giratoire irrotationnel possédant la méme constante de gauchissement C 
et le méme nombre de Mac# local M; au rayon intérieur. En combinant les 
relations (la) et (12a) avec les équations fondamentales: 


LS ee tea Beare (17) 


2 
et 
p THEY) p* T#-klk-1) (18) 


on obtient la relation correspondant a (16): 


M? 

Zee we) =e 2 2 
y ae 1 4 (1+ C2) y 
Ce ye Z : Pees ys 


M*Xy) = (19) 


La figure 2 reproduit les graphiques des fonctions de Macu définies par les 
équations (16) (écoulement supersonique devant l’onde de choc), (10) et (16) 
(écoulement subsonique derriére l’onde de choc) et (19) (écoulement irrota- 
tionnel) pour «; = 30° et M, = 1,363, ce qui correspond a P; = 2 avec k = 1,48). 
Bien que cette comparaison permette d’ores et déja d’affirmer que |’écoulement 
transsonique vérifie la condition de stabilité radiale des deux cétés de l’onde 
de choc, nous n’établirons définitivement cette propriété qu’aprés avoir étudié 
la distribution radiale d’énergie correspondant a (16). 

A cet effet, nous supposerons que l’écoulement giratoire défini par (16) et 
(la) est engendré depuis le repos [état d’arrét initial uniforme (6%, T*)] par un 
ou plusieurs étages préparatoires qui effectuent sur le fluide un travail bien 
déterminé et fonction du rayon, que nous allons calculer. L’état d’arrét initial 
et l’état d’arrét (p*, T*) du fluide devant l’onde de choc sont liés par l’équation 
de l’énergie: 


A 
Diss Tips _ A(y) , (20) 


1) Débit rapporté a lunité de section frontale de l’étage. 
*) Nous verrons plus loin qu’il y a lieu d’envisager encore d’autres limites au rapport caracté- 
ristique v. 


3) Cet exemple numérique ne représente pas nécessairement un cas d@ application pratique. 
Il mettra surtout en relief l’importance des influences radiales étudiées ici. 


5 
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ot H(y) est le travail effectif (kgm/kg) en question. L’évolution thermodyna- 


_ mique du fluide peut étre décrite par la relation polytropique: 


p* Tx ml (m1) _ px Lepr (21) 
qui lie les états d’arrét initial et devant l’onde de choc. De plus: 


m k 
mene lps (22) 


S’il est indiqué d’admettre pour & une valeur moyenne constante pour tout 


écoulement, il est nécessaire de tenir compte de la variation du rendement 


“N 
So 


x 


M,M 
sw 
% 


S 
S 


Nombres de Mach Jocau x 
S 


SS 
S 


4 
Ordonnée radiale y= Tr 


Figure 2 
Nombres de Macu locaux des écoulements (16), (10) et (19). 


polytropique 7, sur le rayon. Nous avons montré ailleurs l’importance de ce 
fait pour les compresseurs axiaux subsoniques [6]. Dans le cas qui nous occupe, 
un calcul de H(y) avec 7, constant conduirait 4 une détérioration inadmissible 
de l’onde de choc, consécutive a celle du profil des vitesses 4 son amont. Nous 
verrons d’ailleurs que l’admission d’une fonction 7,(y) ne complique nullement 
le calcul. Les états statique et d’arrét du fluide devant l’onde de choc sont liés 
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par les relations (17) et (18), qui permettent d’écrire: 
LEED k—1 4 o)ke-1) 
= [i+w], 


de sorte que (21) peut se mettre sous la forme 


k k/(k—1) 


p [r+ 2e+ me = 95 [1 el (21a) 


ou nous avons posé: 


Logarithmons, puis dérivons (21a) par rapport a y et introduisons la condition 
d’équilibre radial (12a) devant l’onde de choc: 


nae 20 A g(t = AG 


P (23) 


SS * [ip log (1 + sx)] 


Le membre de gauche de cette équation différentielle est une fonction connue 
de y, donnée par (16) et (1a). Les variables étant séparées, (23) s’intégre par 
simple quadrature. On obtient ainsi: 


eet 2(y v 24) 
1+ M (v) : ‘ Z ( 
+ log | - a (pk y / M*() ae a(&) dé 
i —{— M} , § 
et, finalement : 
Ral 1/np(y) 
AT 2 
Ter HY) [1 ain mle hit 2 nite 
te CEL Te a} k 


foe M2 
: " (25) 
k=1 " M*(E) cos? «(&) 


xX EXP ’ 
“y(¥) i 


ae, 


E 
S 


qui montre l’importance de la fonction 7,(y) et la simplicité avec laquelle elle 
s’introduit dans le calcul. Il est 4 peine besoin de remarquer que la formule (25) 
est générale et qu’elle s’applique 4 tout écoulement giratoire compressible. 
Nous donnons sous forme de tableau quelques détails numériques du calcul 
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du rapport H(y)/H¢ pour l’exemple considéré plus haut, que nous spécialisons 
a cet effet pour y,=2 (v=0,5), T*# = 288°K, c, = 0,240, c’est-a-dire 
H¢ = 29530 m et H; = 6000 m. 


y | 1 | 1,25 | 15 | 1,75 | Q 
5 COS?a 
co (391, . \ 0793. | -0:519,. | 20.364. | 0,266 
na 0 0,259 | 0,413 | 0,527 | 0,598 
Np(¥) 0,800 | 0,895 | 0,900 | 0,885 | 0,710 
A; ]"p,/nv\v) 
ii: ral 1,2032 | 1,1798 | 1,1787 | 1,1820 | 1,2316 
0 
ag) Unpls) 
1 + —~— M? 
2 (y) : 1 1 1 1 
p= 
ees a 1,0605 | 1,0870 | 1,1033 | 1,1417 
Hiya 0,203 0,249 
ee 0,303 | 0,359 | 0,511 
0 


Cet exemple numérique montre bien l’importance des influences radiales. 
Nous en tirerons deux conclusions: ; 

a) L’écoulement giratoire transsonique est stable. En effet, le travail H(y) 
croissant avec le rayon, il en est de méme de la température d’arrét au droit 
de l’onde de choc. Mais le nombre de Macu local ne diminuant pas du rayon 
intérieur au rayon extérieur dans la méme mesure que pour l’écoulement 
irrotationnel, le drall de l’écoulement supersonique (16) augmentera avec le 
rayon: cet écoulement est donc stable. A plus forte raison, l’écoulement sub- 
sonique en aval de l’onde de choc est stable, car son nombre de Macu local 
croit avec le rayon. 

b) Supposons que la distribution radiale H(y) du travail préparatoire soit 
réalisée par un seul aubage moteur. Si H était constant, l’écoulement en amont 
de cet aubage serait irrotationnel, en vertu de son énergie d’arrét constante et 
de la fonction de gauchissement potentielle (1a) 4 la sortie de l’aubage. H n’étant 
pas constant, le drall en amont de l’aubage moteur préparatoire ne le sera pas 
non plus. Mais il ne pourra décroitre du rayon intérieur au rayon extérieur, 
car l’écoulement giratoire correspondant serait instable. Dés lors, il est facile 
de voir que le drall en question doit étre dirigé en sens inverse de la rotation 
de la machine, en vertu de la relation fondamentale d’EULER 


woAd 


og 
5 


H= 


ZAMP V/10 
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et de la fonction croissante H(y). Dans ces conditions, seul un aubage moteur 
préparatoire supersonique entre pratiquement en considération pour engendrer 
l’écoulement supersonique envisagé, d’une part a cause du contre-drall et de 
la vitesse périphérique élevée de l’aubage moteur, d’autre part a cause de la 
variation considérable du travail H(y) sur le rayon, qui engendrerait dans un 
aubage moteur subsonique une distribution radiale inadmissible des charges 
aérodynamiques. 

Dans l’exemple que nous venons de calculer, le travail de l’aubage moteur 
préparatoire est 2,5 fois plus grand au rayon extérieur qu’au rayon intérieur. 
C’est évidemment une difficulté qui se manifestera dans le dessin de l’aubage 
et qui fera apparaitre une limite aérodynamique au rapport caractéristique ». 
Il y a deux moyens efficaces de réduire la différence H, — H;: 

1° Réduire la hauteur relative des ailettes de l’aubage directeur trans- 
sonique, ce qui diminuera sa section annulaire et son débit spécifique, mais pas 
nécessairement celui du compresseur, qui dépend avant tout du rapport » du 
premier aubage et de la vitesse d’entrée du fluide. 

2° Augmenter l’angle «, a l’entrée de l’aubage directeur transsonique, ce qui 
réduira le drall moyen devant l’onde de choc et l’importance des influences 
radiales, mais également le travail de l’aubage moteur préparatoire et le taux 
de la compression transsonique. 

Dans la plupart des cas, il s’agira donc de chercher une solution de compro- 
mis. Nous situerons la question en traitant un deuxiéme cas limite: celui ou 
l’écoulement supersonique qui aborde l’onde de choc frontale est purement 
axial. Les influences radiales disparaissent et l’on voit immédiatement qu’il 
n’existe aucune raison d’engendrer une distribution radiale non uniforme d’éner- 
gie mécanique. L’écoulement supersonique en amont de |’onde de choc est donc 
irrotationnel, abstraction faite de l’influence de la répartition 7,(y), et il ne 
peut étre engendré que par un aubage moteur supersonique qui redresse l’écoule- 
ment a contre-drall produit 4 son amont par le premier aubage directeur de la 
machine. En principe, l’aubage moteur préparatoire est du type a action; 
pratiquement, il sera avantageux de prévoir une légére réaction, qu’on réalisera 
en adaptant. en conséquence la hauteur des ailettes 4 la sortie de l’aubage 
moteur. On réalisera ainsi un étage transsonique a degré de réaction négatif, 
dont les possibilités sont rapidement évaluées: 

Admettons 


Vy = O85 et C1, = 0,9 ay, = 290 m/s. 
Choisissons 
%, = 30°, cest-a-dire ¢,,, = 251 et lmi, = 145, 
partant 


%, = 49,1°, c’est-a-dire cj, = 125,5 et Cmig = 145. 


ee a ee 
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Choisissons encore 
u; = 200, ce qui donne pi, =u lhe Sc. 
et 
U%q_= 400, ce quidonne f,, = 15,4°. 
Les valeurs de f,, et 6, peuvent étre pratiquement considérées comme limites 
inférieures. On obtient ainsi 
@,,=475m/s et w,, = 546 m/s. 


L’écoulement relatif dans l’aubage moteur est supersonique. Le travail moteur 


vaut: 
_ 4 Ac 200.251 
H = 6.81 5120 m, 


c’est-a-dire 
H =42 kcal/kg .denc AT(H)= 50°C, 


Quant a la vitesse de l’écoulement qui aborde l’onde de choc, elle serait ap- 
proximativement de 560 m/s. On aurait donc M ~ 1,8 et, d’aprés (7): 


ee De 3 6) 


Parlor 2,4 


ce qui correspondrait a un rendement de compression 7, ~ 0,835. En admet- 
tant un rendement adiabatique total de |’étage de 0,75, le taux de compression 
total serait 

i). 31,0, of © =i 955 (CINA 0). 


On voit le désavantage du procédé: Le travail moteur est limité par la con- 
dition d’entrée axiale dans l’aubage directeur transsonique et par les angles 
d’entrée relatifs dans l’aubage moteur. L’expansion nécessaire a la production 
de l’écoulement supersonique axial est compensée par un choc intense de 
rendement médiocre. Il n’en résulte qu’un taux de compression relativement 
faible pour |’étage entier. On pourrait se contenter d’engendrer un écoulement 
relatif sonique ou légérement supersonique devant l’aubage moteur. Mais, pour 
un débit donné, cela reviendrait 4 réduire le travail moteur réalisable et par 
conséquent, le taux de compression de 1’étage. 

Il est donc nécessaire d’utiliser un écoulement transsonique giratoire pour 
réaliser dans des conditions favorables un taux de compression P* voisin de 2, 
valeur qui, 4 notre avis, correspond a ce qu’on devrait attendre d’un étage de 
compression par onde de choc frontale. 
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Nous pensons que les deux exemples esquissés situent bien la question de la 
compression transsonique dans l’aubage directeur et montrent suffisamment 
l’importance des formules (16) et (25). 


Remarques 


Dans le calcul de la distribution radiale H(y), nous avons admis que l’équa- 
tion polytropique (21) lie les états d’arrét initial et devant l’onde de choc. 
Il en est résulté des calculs simples. Mais il va de soi qu’on aurait tout aussi 
bien pu détailler l’évolution thermodynamique du fluide en admettant un 
rendement polytropique 7,, pour l’expansion dans le premier aubage directeur 
et un rendement polytropique 7,, pour la compression (états d’arrét!) dans 
l’aubage moteur. Sans étre inextricables, les calculs en seraient cependant 
devenus fastidieux, raison pour laquelle nous n’avons pas jugé opportun de 
les reproduire ici. Remarquons seulement que notre rendement polytropique 7, 
est ainsi fonction des rendements 7,, et 7,, et qu’il est évident que son admis- 
sion devra si possible étre étayée par l’expérience. 

Dés le début des calculs, nous avons introduit des variables sans dimension, 
en particulier des nombres de MAcH locaux. Ici également, ce choix a été dicté 
par les simplifications qui en résultent dans le développement des formules, 
plutot que par souci d’esthétique. On reviendra aisément aux grandeurs di- 
mensionnelles a l’aide de la formule (20a) et des relations fondamentales. Nous 
n’€crirons ici que les expressions de la vitesse locale du son en amont et en 
aval de l’onde de choc: 


ile ue Bee 
a ay en et (26) 
eee 
2 
et 
. i, k-1 
7S. k M? — —~— 
an) ee (k—igH 2 
i= Vr am |/ p+ PSR hearer (27) 


ou ag est la vitesse locale du son & la température d’arrét initiale Tx. 

Enfin, nous n’avons donné aucune indication sur le calcul de la distribution 
radiale des vitesses devant l’aubage moteur. Leurs composantes périphériques 
sont déterminées par la condition (1) et la fonction H(y) en vertu de l’équation 
fondamentale d’EuLER, mais il y aurait lieu de tenir compte des déplacements 
radiaux des lignes de courant a la traversée de l’aubage moteur, qui risquent 
d’étre considérables a cause de la distribution radiale non uniforme de énergie 
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du fluide devant l’onde de choc. Nous avons abordé ce probléme dans 1’étude 
[6] déja citée. Nous n’y reviendrons donc plus. 


III. Conclusions 


Dans le cadre d’un exposé destiné a établir quelques principes généraux, 
nous n’avons discuté que des problémes directement accessibles 4 l’étude ana- 
lytique. Les conditions d’écoulement étudiées sont nécessaires A la réalisation 
d’un compresseur axial supersonique 4 ondes de choc frontales possédant un 
rendement de compression élevé. Elles ne sont nullement suffisantes, car elles 
n’assurent pas la stabilité des ondes de choc. Cette condition importante n’est 
satisfaite que si chaque élément de flux annulaire évolue, en aval du choc, 
vers une section de sortie de l’aubage bien déterminée, fonction des conditions 
d’écoulement au droit de l’onde de choc, partant du rayon. C’est un probléme 
dont la solution devra étre recherchée pour chaque cas particulier et qui sera 
d’autant plus facile a résoudre que les paramétres de |’écoulement varieront 
moins avec le rayon d’aubage, c’est-a-dire que l’intensité de rotation du fluide 
sera plus faible. De toute facon, il en résultera une servitude supplémentaire 
dans le dessin d’aubage et, cas échéant, une nouvelle limite aéro-thermodyna- 
mique au rapport caractéristique v= ~7,/r;, qui dépendra, entre autres fac- 
teurs, des coefficients de charge aérodynamique admissibles de l’aubage. II 
conviendrait également de comparer les possibilités du compresseur a ondes 
de choc frontales dans l’aubage directeur et celles des autres types de com- 
presseurs axiaux supersoniques. On trouvera quelques prévisions a ce sujet 
dans les études [1] et [2]. Nous nous bornerons a rappeler ici une difficulté 
inhérente 4 la compression transsonique dans l’aubage directeur, a savoir les 
frottements élevés du fluide dans l’aubage moteur, dus a l'utilisation d’un 
écoulement relatif supersonique. 

On a insisté récemment sur les difficultés rencontrées dans la stabilisation 
de l’onde de choc dans l’aubage moteur [7]. Ces difficultés auraient méme 
motivé l’abandon de ce mode de compression au profit de celui étudié ici. I] est 
-malheureusement difficile de savoir si l’on a réalisé dans ce cas les conditions 
nécessaires a l’existence d’un écoulement relatif transsonique stable: dans 
maints pays, |’étude des compresseurs axiaux supersoniques a un coté «Défense 
nationale». On en est ainsi réduit a des conjectures. 

On a également insisté sur l’épaississement des couches limites a la traversée 
de l’onde de choc et sur les difficultés qui en résultent. La succion ou le soufflage 
éventuels de ces couches limites est évidemment plus facile a réaliser dans 
l’aubage directeur que dans l’aubage moteur. Sans minimiser l’importance de 
ce probléme, nous pensons toutefois qu’en tenant compte de la distribution 
radiale non uniforme du rendement 7, dans le dessin d’aubage, on réduira 
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dans une large mesure une difficulté qui semble aujourd’hui prohibitive. Quoi 
qu’il en soit, l’abandon de la compression transsonique dans l’aubage moteur 
ne peut se justifier que par V'impossibilite de réaliser des conditions d’écoule- 
ment adéquates, c’est-a-dire analogues a celles que nous venons d’étudier pour 
le cas de l’aubage directeur transsonique. Comme nous l’avons dit plus haut, 
c’est une question sur laquelle nous nous proposons de revenir plus tard. 


Je tiens & remercier Monsieur F. Rotu, Directeur de la Société Ardag S.A., 
pour son autorisation de publier cette étude. Qu’il me soit également permis 
d’exprimer mes remerciements 4 Monsieur le professeur J. ACKERET, qui a bien 
voulu s’intéresser a ce travail. 
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Summary 


In an attempt to study the three-dimensional flow in a supersonic axial 
compressor, the normal-shock-in-stator case has been considered. It is shown that 
the kinematic and dynamic conditions which are necessary to produce a normal 
shock wave and a radially stable transsonic flow are very restrictive: 

(i) The radial distribution of the flow angles in front of the shock wave is found 
to be that of the classical irrotational giratory flows. 

(ii) The radial equilibrium condition of the flow on both sides of the shock wave 
requires a well determined radial distribution of the MacH numbers of the 
velocities approaching the shock, which results in a non-uniform radial energy 
distribution of the transsonic flow. 

Furthermore, it is necessary to take in account the radial variations of the 
polytropic stage efficiency 7, in order to induce the prescribed energy distribution 
through one or several shock-free rotor-blade rows. A numerical example shows 
the importance of the radial effects. Finally, some critical comments are made 
on the shock-in-rotor case. 


(Regu: le 6 mai 1953.) 
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Note on the Boundary Layer on a Rotating Sphere 


By Swami Daya. Nicam, Kharagpur, India?) 


1. Introduction 


HowartTH?) has investigated the problem of the flow engendered by a sphere 
rotating uniformly about a diameter in otherwise undisturbed fluid. He considers 
the sphere to be made up of two hemispheres joined smoothly at the equator. 
The flow being symmetrical about the equatorial plane, the boundary layers 
originate at the poles and develop towards the equator where they impinge on 
each other. He has obtained a solution of the boundary-layer equations in the 
form of power series and has shown that the flow at the poles approximates to the 
rotating disk solution of voN KARMAN (1921). Due to lack of rapid convergence 
this solution cannot be used in the vicinity of the equator. Even the approximate 
solution of the same problem by the Karman-Momentum Integral Method does 
not give satisfactory results near the equator, because some of the assumptions 
involved break down. The solutions are defective in so far as they do not give 
any indication of outflow of the liquid near the equator. HowartH attributes 
the cause of the failure of his solutions near the equator to the boundary-layer 
equations, which, he says, must fail to represent the region of interaction between 
the two impinging layers on account of the parabolic character of the equations. 

We, however, do not agree with this explanation, because the assumptions 
of the boundary-layer theory do not break down at the equator, and think that 
the equations are quite sufficient to describe the interaction at the equator. 

In this note we have shown that it is possible to construct a solution in power 
series form which fulfils all the physical requirements of the problem. It is difficult 
to say anything definite about the convergence of the solutions near the equator. 
But the solutions have a definite advantage over those of Howartu, because 
apart from the inflow at the poles they give a clear indication of the outflow 


near the equator. 


2. The Boundary-Layer Equations 


In the following treatment y, 0 and @ represent the spherical polar coordinates 


_ with y measured radially outwards from the centre of the sphere, 6 measured 


from the axis of rotation and » the azimuth. w, u, v represent the velocities in the 
directions of 7, 6 and @ increasing respectively. The equations of motion in 
spherical polar coordinates with terms of azimuthal variation omitted, are 


Ow Ow r u dw yu? + y2 
2 eS = Ww a — . — = “ 
OF Oy 00 y 
0 . : 2 2 0 2 té A) 
1 2) Hie, oe gt uU oe] 
= 7 hein enya ae aa 


1) Department of Applied Mathematics, Indian Institute of Technology. 
2) L. Howarrn, Phil. Mag [7] 42, 1808-15 (1951). 


Ou Ou u Ou uw v2 
de : sat =e coté 
a Gt Oey | 
1 Op Dine OL take eed 
=— ay ont? [PM oe ae — Frain t (A) 
ov Ov Cee) vw Uv a [ Cae rusty | 
GT” in kak a aaa ial A | 
1 ) 1 0 : 
ee TERE Bee Ty == Ol, 
y2 ” Oy (y? w) + ysin@ 00 eee) 


To obtain the boundary-layer equations, we take w ~ 0(6), uw ~ 0(1), v ~ O(1), 
and (0/0y) ~ 0(6-}). The boundary-layer equations for a sphere of large radius are?) 


u Ou Ou v2  02u 


Keepreted ee eee Gt ean ee 1 
rear) +w i Pao a dy? (1) 
u Ov ov uv 0*v 

~ eecen = 2 
reeks was cord =» 5, (2) 
1 ou Ow u 

é3 a) — . 5 
ee ge cot = 0 (3) 


It is clear that the assumptions of the boundary-layer theory do not break down 
at the equator. 

As remarked earlier that the flow near the poles approximates to that due 
to a rotating disc, it is interesting to note that the flow near the equator approxi- 
mates to the flow due to a cylinder rotating about its axis and the liquid stream- 
ing along the axis from both sides. In the boundary-layer equations, if we neglect 
the terms containing cot @ (cot 0 > 0 as 0 > 2/2), we get 


u Oe tel py Olas ago 3O Oth 4 
a 00 Oy 4 Oy?’ (4) 
u ~ Ov ia Ov ne: 02u 5 
a 00°” dy” Ody? (5) 
1 Ou Ow 

@ Mag Gakgg aoe (9) 


These equations are identical with the equations derived by HowartTH in the 
paper cited in the introduction. The solutions of these equations involve an out- 
flow of the liquid. Hence the equations (1), (2), and (3) seem to be quite sufficient 
to describe the flow near the equator. 

For the problem of the rotating sphere we take 


u=aQ2cos6 [sind F, + sin?6 F, + sin'é K+ .--], (7) 
v =a Q{(sin@ G, + sin’6@ G, + sin'@G, + ---], (8) 
5, y 2\1/2 : ase Aas - 
w= (=) (3 cos?@ — 1) [H, + sin?@ H, + sint@ H, + ---), (9) 
Q\1/2 
ay was (~) (10) 


1) L. Howarrn, Phil. Mag [7] 42, 1808-15 (1951). 


a 
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where a is the radius of the sphere and 2 the angular velocity. We get the 
following equations for the first nine functions: 


oer ee a Oe ee (11) 
2G MTP a mon. (12) 
2Fe tHe 0., (13) 

’ y 3 
eee Ret Pe oO, G2 HE Bh (14) 
22, Gt 4 Fy Get Hy Ge Hy G4 2G, Hy Gh= GP, (15) 
4F,4+ H,=0, (16) 


6FRR—6F7,i4+3F}4+H,Fi+FiH,+ FH, 
3 / f mu” 
— > (Hy Fy + Fi H,) — G}- 26,6, = Ff, 


Mirahiet, Pat 4G 4 GF 2G, e+ Hyp? 


3 , , s , W 18 
ey (H, G, + G, H;) + G3; H; + G, H, = Gj, 


6F,+ F,+ Hj=0. (19) 
The boundary conditions are 
FiG=i,=—Fo= 0; G,=1, Gz=G;=—0; H,=H,=—H,= 0atZ=0. 
F,>0,F,>0,R>9; G,>0, G,>0, G, > as Z > o~. 


OL a Ne ee eee ee 


§ The equations (11) to (19) have been solved approximately by the method 
_ of PoHLHauseEeN. The functions F and G have been taken as polynomials in s 
3 where Z = s 6 and the functions H have been calculated from the equation of 
4 continuity. 

; We take 

4 1 1 

BP F=as(l—s*(1+2s)—-Fst(t—s)*, Gy= = (2+ 5) (1-3), 

‘2 

: if 

Pe s(t —“s)*{1 + 2 sy s* (1 — 8)", (COR (Ee See 

fy 

y 1 2 2 e 2 

: Poa ae (b— s)*(1 + 28) — | s2(1—s)*, Gaees (i —.:5)2., 

se 


The functions have been chosen to satisfy the following conditions: 


F,=F,=F,= 0; G,=1, G,=6,= 0; Hj, = H,=H;=0 on s=0o0rZ=0. 


Aero a 


F,= F, =F, = 0; Gy=G,=G,=0 on s=1 or Z= 6, 


6 being the thickness of the boundary layer. In order that the flow b> continuous 
and smooth on the outer edge of the boundary layer, we should have, 


F{(1) = Fg(1) = F5(1) = 0 and G4(1) = G5(1) = G5(1) = 0, 


where a dash denotes differentiation with respect to s. 


ee AE ees ee hh ee aE 


— 
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Integrating (11) and (12) between the limits 0 and 6 using (13) we obtain 


ul 1 1 
; 3 
3 Fids— [Gids=-ad*, 4 | F, G, ds = 
0 0 0 


Figure 1 
The flow functions F, G, and H. 


Whence, after substituting the values of F, and G,, we get a = 1.037 and 
6 = 2.5587. The equation of continuity (13) gives 


: ae E 1 Se s® 
H, = —5.1174 [1-037 (0:5 s?— 0-75 s4 + 0.4 5) — — st + | 


Similarly, for the remaining equations we find 


b= 0.4715; c = 05997; d= 0.3440; e = 0.3478 ‘ 


53 54 5) 
H, = —10.2348 [o-4715 (0.5.57 — 0:75 ¢8 se 0.4 58) | ? 


H 1 s3 s4 $§ 1 
5 = —15-3522 | 0.3440 (0-5 s? — 0.75 s4 + 0.4 55) — te ee 78 Set Fe 
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The numerical values of the functions are tabulated below: 


Table 1 


s | O-1 | 0:2 | 0-3 | 0-4 | 0:5 | 0-6 | O7 | 0-8 | 0-9 | 1-0 


0-0976 | 0-1730 | 0-2218 | 0-2400 | 0-2280 | 0-1902 | 0-1347 | 0-0735 |.0-0220 0 | 
0-8505 | 0:7040 | 0-5635 | 0-4320 | 0-3125 | 0-2080 | 0-1215 | 0-0560 | 0-0145 0) 
0-0225 | 0-1010 | 0-1974 | 0-3173 | 0-4383 | 0-5461 | 0-6299 | 0-6802 | 0-7063 | 0-7108 
0-0417 | 0-0717 | 0-0888 | 0-0934 | 0-0866 | 0-0708 | 0-0492 | 0.0264 | 0.0078 0 
0-0486 | 0-0768 | 0-0882 | 0-0864 | 0-0750 | 0-0576 | 0-0378 | 0-0192 | 0-0054 0 
0:0226 | 0-0868 | 0-1645 | 0-2590 | 0-3523 | 0-4329 | 0-4951 | 0-5337 | 0.5496 | 0-5537 
0-0293 | 0-0488 | 0-0588 | 0-0604 | 0-0547 | 0-0493 | 0-0299 | 0-0158 | 0-0046 0 
0-0282 | 0-0445 | 0-0511 | 0-0501 | 0-0483 | 0-0334 | 0-0219 | 0-0111 | 0-0031 0) 
0-0300 | 0-1130 | 0-2098 | 0.3260 | 0.4392 | 0-5347 | 0-6080 | 0-6523 | 0.6700 | 0-6750 


SO OO 


i) 


3. Discussion of the Results 


The expressions (7), (8), and (9) satisfy all the conditions on the surface of 
the sphere. Both w and v tend to zero on the outer edge of the boundary layer; 
u vanishes on the equatorial plane 0 = z/2 inside the boundary layer; w is nega- 
tive so long as 3 cos?@ — 1> 0, it vanishes for 3 cos?@ — 1 = 0, i.e. 6 = 54° 45’ 
in the upper hemisphere and 9 = 125° 15’ in the lower hemisphere, and is positive 
when 3cos?@ — 1< 0. This shows that there is a radial inflow in the upper 
hemisphere for 0 < 0 < 54° 45’; this inflow disappears on the cone 3 cos?@ = 1 
and changes into an outflow for 54° 45’< 90 < 2/2. Similar conclusions can be 
drawn for flow in the boundary layer on the lower hemisphere also. Nothing 
definite can be said about the convergence of the series solutions near the equator. 


Zusammenfassung 


Die gleichmassige Str6mung in der Grenzschicht bei einer gleichmAssig rotie- 
renden Kugel wird behandelt. Die Grenzschicht entsteht an den Polen und ent- 
wickelt sich in beiden Hemispharen nach dem Aquator zu. An den Polen verhalt 
sich die Kugel wie eine rotierende Scheibe, und Fliissigkeit str6mt hier ein. Am 
Aquator verhalt sich die Kugel wie ein rotierender Zylinder, und die Flissigkeit 
str6mt in dessen Bereich aus. 


(Received: May 13, 1953.) 
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Boundary Layer on Rotating Spheroids 


By Buaskar SADASHIV FaApNnis, Kharagpur, India?) 


1. Introduction 


HowartH?) has attempted the solution of the boundary-layer equations of 
the flow engendered by a sphere rotating uniformely about a diameter in an 
otherwise undisturbed fluid in the form of an infinite series. He finds that bound- 
ary layer originates at the poles and develops towards the equator. His solution 
approximates to the rotating-disc solution of von KARMAN’) at the poles but does 
not give an outflow near the equatorial region. HowartTH attributes the cause of 
the failure of his solution near the equator to the boundary-layer equations, which 
according to him must fail to represent the interaction between the two imping- 
ing layers at the equator on account of their parabolic nature. 

Nicam‘) has shown that it is possible to construct a solution in power-series 
form, which fulfils all the physical requirements of the problem. He thinks that 
the boundary-layer equations are quite sufficient to depict the interaction be- 
tween the two impinging layers because the assumptions of the boundary-layer 
theory do not break down at the equator. Nothing definite, however, can be 
said about the convergence of these series’ solutions near the equator. 

In the present paper the flow generated within the boundary layer by the 
rotation of prolate and oblate spheroids have been separately considered and the 
results compared with that of a sphere having the same equatorial circle. 


2. Boundary-Layer Equations 


In what follows «, B, y represent the orthogonal curvilinear coordinates, and 
“u, VU, W, are the components of velocity in the direction of a, B, y respectively. 
The equations of motion in orthogonal curvilinear coordinates with terms of 
azimuthal variation B omitted are 


w {0 ) vo (0 \ 1 Op 
een Wye, Vee plies rt eB as 
Pe ae al “hs wis aoe oh, Oa 


= oot aloe (“3 hs — _ 9 (hgh, a il 0 (/ur+tv2+ w2 


i (2.1) 
w U 0 
Hescpiae he v) 4: speeee iis 
RE Ronee | 
CGA) 
Pal Ee (= hs -) 3 0 (=: uid dv 
hy Ba By On Naas Oe. Fy hea Ole al =| ‘ 
u 0 0 v 0 il Op 
- -(h, u) — a (he 2 + ——___ . - 1, v) = —— - - 
hy hs Ee wre) 0a (is w)| " Noh, oy (2 ) oh, Oy 
LY oh [5 (4 hs 3) " 0 (> ho Ow 4 il 0 (u®+ Ee 
iy highs 0x \ hy 00 Oy\ hs arti hs Oy ({ =) : 
(2. 3) 
1) Department of Mathematics, Indian Institute of Technology. 
*) L. Howarrn, Note on the Boundary Layer on a Rotating Sphere, Phil. Mag. [7], 42, 1808-15 
(1951). 
3) T. von KArmAn, ZAMM. 1, 244-47 (1921). 
AESsD: pee Note on the Boundrs Layer on a Rotating Sphere, ZAMP 5, 151 (1954). 
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and the equation of continuity is 


il ie eR 


i 1 0 i) 
Tc tirTs Yow (ia Bs) + Ze (hah w)| = 0, (2.4) 
4 ‘where h,,/2,h, are the length elements in the directions «, 8, y. Following 
PRANDTL’S argument, we take u~ 0(6). v~ 0(1), w~ 0(1), 0/da~ 0(6-1), and 
assuming’) that h2, 3, and their derivatives are of ~ 0(1) the equations (2.1), 
(2.2), (2.3), and (2.4) reduce to 


q 


1 0p 
acta 0(1) or pro(d), (2. 5) 
4 w ov, u_ ov wu Ohs yO 2.6 
hy On Ty Ot . lights . Oy he- O03” (9) 
wow, w ay vt (Oh, - v  O*W (2.7) 
hy- 0% hy OY © hah, Oy HE 008” 
1 Ou 1 ow w 0 
pete oh by ae io 


3. Prolate Spheroid 
For a prolate spheroid we introduce a system of coordinates defined as follows: 
Z+tr=kcosh(a+tzy), gp=68, (3.1) 


where +, z, m are the cylindrical coordinates and «, y are the elliptic coordinates 
_ in the meridian plane, and f is the azimuth. After substituting 


hy = hs = k Vcosh?« — cos*y, h, =k sinha siny 


and writing « = «) for large spheroids, equations (2.6), (2.7), and (2.'8) become 


Ov ov v 07v 
“u— + w—+wucoty = a (3. 2) 
0 oy eg Vcosh?a)— cos?y 9? 
Ow Ow v 0?w 
u—— + w—— — v* coty = — ———— + (353) 
0a oy ane Vcosh?a) — cos?y 9%" 
Ou Ow cosy siny 
os ea igs ean ot) "9% 3.4 
Om + Oy aoe [coty e cosh?a, — | Weg 
For the components of velocity”) we take 
1/2 
7 (7°) (3 cos*y — 1) [H, + sin’y H, + ---], (3. 5) 
v=hkQsinha, siny [G, + sin?y Gg +---], (3. 6) 
w= kQsinha, cosy siny [F, + sin’y Fy + ---], (3.7) 
1/2 
a (=) : Rk sinha, (% — &) . (3. 8) 


- 


1) L. Howartu, Phil. Mag. [7], 42, 239-43 (1951). 
2) S. D. Nica, Note on the Boundary Layer on a Rotating Sphere, ZAMP 5, 151 (1954). 
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For the first six functions we get the following six equations 


FP? G@4+.Fi = Fis (3.9) 
2 Gi Hy Gla Gl (3.10) 
a» in pe ye 
, , gs! aa E ah apheed Fi 
47, F,+ Fi H,+ H,F,—26,G6,—2Fj— > HW, Fy = eS beeen 
(3.12) 
/ / 3 / nN 1 Gi 
ZF Gut: 4 FG, Hy Gy cb By Gem 20 Gy ee, Go ae nha 
(3.23) 
4F, i ee = (3.14) 
Sa Feri “a, 


The boundary conditions for these functions are 
yeh eens heen) Co ne Cree es elie Wel (Vache Bee nl ere Sie alte (0 
G,>0, G,>0atz> o. 


The equations (3.9) to (3.14) have been integrated numerically by the method of 
POHLHAUSEN. The functions F and G have been taken as polynomials in s where 
z= s 0, 6 being the thickness of the boundary layer. The functions H have been 
calculated from the equations of continuity. 


R= as (1—s* +25) — 5 hae ce (3.15) 
G=52+9 1-9), (3.16) 

Fie bail org Peis) eae eh aa (3.17) 
Go e's. (1 — 5) (3.18) 


These have been chosen to satisfy the boundary conditions 

h= = 0,°G,= 1 G,= 0, A,—H,= 0 
on s equal to 0 or z= 0 

G=G,=0, F.=R=0 on s=1-orz= 6. 


In order that flow be continuous and smooth on the outer edge of the boundary 
layer we should have 


Fi(1) = Fi(1) =.0 and G{(1) = Gi(1) = 0 


where a dash denotes differentiation with respect to s. Integrating (3.9) and 
(3.10) between the limits 0 to 5 and using (3.11), we obtain 


1 1 
3 | Fras | Gjas——ao-, + [ FG, as = 4 é-2, 
0 0 


where after substituting the values of F,, G, we get a = 1-037, and 6 = 2.5587. 
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The equation of continuity (3.11) then gives 
ee i) 


Taking sinh?a, = 10 we similarly obtain from the equations (3.12), (3.14), (3.17), 
and (3.18) 6 = 0.4730, c = 0.5299, and 


z 1 Se Se 
H, = —5.1174 [1-037 (0-5 s? — 0-75 s# + 0-4 s®) — - gs? 4+ —_ — ; 


‘ 3 4 5 
eH, = —10-2348 [0.4730 (0-5 s? — 0-75 4 + 0.4 58) — < ~ < = 7 + 0127935 H,. 


! 


It is interesting to note that the equations (3.9), (3.10), and (3.11) are the same 
as in the rotating-disc problem of von KARMAN. 
4. Oblate Spheroid 

For the oblate spheroid we introduce a system of coordinates defined by 
zg+ivy=ksinh(«+ty), o—B,; (4.1) 

where again /, 2, m are the cylindrical coordinates of a point; « and y are the 
elliptic coordinates in the meridian plane. After substituting in the equations 

(2.6), (2.7), and (2.8), 
h, =h, =k cosh?« — sin?y, h,=kcoshasiny 


10 


Figure 1 
The flow functions F, G, and H. 
a—m— Prolate spheroid; —— sphere; —e~—e— oblate spheroid. 


Pit W == 154453 Pg: %y = 1-8686; O,: % = 1:8686; Og: a = 1.5445, 


"we 
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and writing « = «, for a large spheroid, we have 
Ov ov v 0?u 
“u—+w—-+ wvcoty= . ; (4. 2) 
Oo Oy Crouse cosh? a, — sin?y Oa? 
0?w 
ee oe act me : i“ + = (4. 3) 
Oo Ov k Vcosh?a% — sin? y Oo 
Ou Ow cosy siny 
sei mpecsetgcre On 4.4 
Ou Oy | [eoty cosh?% — Sed a4 
For velocity components?) we take 
1/2 
ain (72) : (3 cos*y — 1) (ta 4+ sin?y Ey -be-), (4.5) 
v=k 2 cosha, siny [G, + sin?y Gg + ---], (4. 6) 


w= k Qcosha, cosy siny [F, + sin?y Fy + ---], 


Z= elk k cosha& (a — a) . (4. 8) 


Figure 2 


The function Gg. 


~a—a— Prolate spheroid; —- — sphere; —e—e~— oblate spheroid. 
Py = G% (15445); Py = o% (18686); O, = % (18686); O, = og (1-5445). 


1) S.D. Nicam, Note on the Boundary Layer on a Rotating Sphere, ZAMP 6, 151 (1954). 
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The equations for F,, G, H, are the same as (3.9), (3.10), and (3.11). For the 
next three functions we get 


‘ 3 1 Fi 
Ree ee Hie 2G. Gi OR? a FT, Re BY 1 
pega fy fig + fy is 1 G3 2 Fi 7 Fi Fy + 2° cosh?a, ’ 
(4. 9) 
ee Grr eae Co Ce 2 mi ce Ow Jy 
3 Gy 1% 3G, + 11, G, 1 sede = s+ 2 * cosh?a, ’ 
{ (4.10) 
AB Ee # = 0 4.11 
at Os coshta, ~ Sia 


The boundary conditions for these functions being the same as in the previous 
case we have assumed for F,, G,, F, G; polynomial expressions similar to those 
used in the case of prolate spheroid. The results of the numerical integration are 
tabulated below. 

Table 1 


| 0-0967 | 0-1730 | 0-2218 | 0-2400 | 0-2280 
0-8505 | 0-7040 | 0-5635 |0-4320 | 0-3125 
0-0255 | 0-1010 | 0-1974 | 0-3173 | 0-4383 


Prolate spheroid: a = 1-8686, 


0-0720 |0-0891 | 0-0938 | 0-0869 
0-0678 | 0-0779 |.0-0763 | 0-0662 
0-0921 | 0-1748 |0-2758 | 0-3755 


Prolate spheroid: a) = 1-5445, 


0-0721 | 0-0893 | 0-0939 | 0-0871 
0-0593 |.0-0681 | 0-0668 | 0-0579 


0-0415 
0-0375 


0-1902 | 0-1347 | 0-0735 | 0-0220 
0-2080 | 0-1215 | 0-0560 | 0-0145 
0-5461 | 0-6299 | 0-6802 | 0-7063 


b = 0-4730, ¢ = 0-5299 


0-0711 | 0-04.94 | 0-0266 | 0-0078 
0-0509. |.0-0334 |.0-0170-| 0-0048 
0-4619 | 0.5287 | 0.5699 | 0.5872 


b = 0:4735, ¢ = 0-4636 


0-0712 | 0-0495 | 0-0266 | 0-0078 
0-0445 | 0-0292 | 0-0148 | 0-0042 


0-7108 


0 
0 
0.5912 


0 
0 


0-0253 |0-0973 | 0-1849 | 0-2920 | 0-3980 |0-4899 | 0-5611 | 0-6047 | 0-6233 | 0-6276 


Oblate spheroid: o, = 1-8686, 6 = 0-4727, c = 0-6637 


0-0937 | 0-0869 | 0-0710 | 0-0494 | 0-0265 | 0-0078 0 
0-0956 | 0-0830 | 0-0637 | 0-0418 | 0-0212 | 0.0060 0 
0-2453 |0-3335 | 0-4095 | 0-4683 | 0-5046 | 0-5194 | 0-5237 


0-0414 |0-0719 | 0-0891 
0-0538 |0-0850 | 0-0976 
0-0261 | 0-0826 |0-1558 


/zcsEEEEnEEE) EEE [ce Se nn a 


Oblate spheroid: % = 15445, b = 0-4726, c = 0.7175 


F, 
| Gs 
— Hy 


0-0710 | 0-0494 | 0-0265 | 0-0078 
0:0689 | 0.0452 | 0-0230 | 0-0065 
0-3888 | 0-4442 | 0-4786 | 0.4924 | 0-4957 


0-0414 |0-0719 |0-0891 | 0-0937 | 0-0868 
0-0581 | 0-0918 | 0-1055 | 0-1033 | 0-0897 
0-0206 | 0-0783 | 0-1483 | 0-2332 | 0-3168 


ZAMP V/11 
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5. Thickness of the Boundary Layer 
Taking air as the fluid, we take 


y= 0-14cm2/s; Q2= 600r/m; % = 1-8686, e= =, @="3 en 


yi1 


since 
ot Qy\1/2 ; 
7 (2:5587) = ey R sinha (a — o%) 5 


thickness of the boundary layer being equal to the distance « — x along the 
normal at any point = 0.00534 cm. 


6. Discussion of the Results 


The expressions for velocities satisfy all the boundary conditions on the sur- 
face of the spheroids. Both u and v tend to zero on the outer edge of the boundary 
layer, w vanishes on the equatorial plane y = 2/2 inside the boundary layer. u is 
negative so long as 3 cos?y — 1> 0, it vanishes for 3 cos?y = 1 or y = 54° 45’ in 
the upper half of the spheroid and for y = 125° 15’ in the lower half of the sphe- 
roid. This shows that there is normal inflow in the upper half for 0 = y S 54° 45’. 
This inflow disappears on the cone 3 cos?y = 1 and changes into an outflow for 
54° 45’< y < a/2. Nothing definite, however, can be said about the convergence 
of the series solutions near the equator. 

In the paper cited in the introduction, NrcAm has shown that the equations 
for the flow in the boundary layer of a cylinder rotating about its axis and the 
liquid streaming past it from both sides can be deduced as a limiting case from 
the boundary-layer equations for the sphere. It is interesting to note that these 
equations can also be deduced as a limiting case from the boundary-layer equa- 
tions for the rotating spheroids. Because the solutions for the cylinder involve 
an outflow in the central plane it may be concluded that the boundary-layer 
equations for the sphere adequately describe the interaction at the equator, and itis 
quite apparent that similar argument also holds in case of rotating spheroids. We 
can in general say that there will be an outflow of liquid in the equatorial plane 
of the rotating sphere and spheroids because near the equator these bodies behave 
like a cylinder rotating about its axis and the liquid streaming in the direction of 
the axis from both sides towards the central plane. The present investigation further 
indicates that the outflow of the liquid near the equatorial region depends 
upon the cylindrical nature of the rotating body near the equatorial plane. 

In the case of spheroids and sphere rotating with the same angular velocity 
and having equal equatorial radii, we find that the value of H, decreases as the 
ratio 7/h increases, y and h having meaning as shown in Figure 3. In Table 2 
values of H; are shown against those of r/h. 


Table 2 
| Prolate spheroids | Sphere | Oblate spheroids 
aX 1-5445 1-8686 1-8686 1-5445 
v{h 0-9129 0-9535 1 1-0506 1:0972 


| — He 0-6276 0-5912 0-5537 0-5237 0:4957 


rey Tye 
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_ The cylindrical surface as shown in Figure 3 is greatest for the prolate spheroid 
_ for which the ratio 7/h is least. For r/h equal to unity we get the sphere and for 
_r/h greater than one we get the oblate spheroid. Therefore the outflow in the 
equatorial plane is greatest for the prolate spheroid and least for the oblate 
_ spheroid while that of the sphere lies between these two extreme values. In the 
_ case of two prolate spheroids defined by 


=) 091297, 1a, = 11-8686, k= 1-05.01, = = 0:9535 


Bo, = 1-5445, k = 1.3488, ; 
L 


i 


0-65 


0-63 


0-45 
a9 092 09% 096 096 1 102 704 106 F108 10 112 
Figure 3 Figure 4.—Graph of Hz against r/h. 


we find that the outflow in the equatorial region is greater for the former. This 
supports the view that the ratio v/f influences the outflow in the equatorial plane. 
Similar conclusions can be drawn for the oblate spheroid also. 

It is of interest to note that azimuthal velocity (Figure 2) at the corresponding 
points in the boundary layer increases as v/h increases. Hence, when two bodies 
of revolution having equal equatorial radii, rotate with equal angular velocities, 
the flatter amongst the two imparts a greater rotational velocity to the fluid 
particles within the boundary layer. 


In conclusion I express my deep sense of gratitude to Dr. S. D. Nicam for his 
kind and unsparing help in the preparation of this paper. 


Zusammenfassung 


Die Grenzschichtgleichungen fiir dreidimensionale Stromung an der Aussen- 
seite eines rotierenden Spharoids werden in konfokalen Koordinaten abgeleitet 
und numerisch mittels der Methode von PoHLHAUSEN fiir verlangerte und ver- 
kiirzte Spharoide gelést. Die Lésungen zeigen einen Zustrom an den Polen und 
einen Abfluss an den Aquator. Der Ausfluss hangt vom Verhaltnis 7/h ab, wobei 
y der Radius des Aquators und h die Hohe des K6rpers in der Rotationsachse ist. 


(Received: May 13, 1953.) 
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Radantrieb und Raketenantrieb 


Von EUGEN SANGER, Paris") 


Man begegnet noch vielfach der Ansicht, der Radantrieb sei dem Raketen- 
antrieb hinsichtlich des 4usseren mechanischen Wirkungsgrades des Antriebs- 
vorganges allgemein tiberlegen. 

Die folgenden, elementaren mechanischen Uberlegungen zeigen, dass dies nur 
fiir den Bereich der uns gelaufigen geringen Fahrgeschwindigkeiten richtig ist, 
fiir die hohen Geschwindigkeiten der Raumfahrt aber keineswegs mehr zutrifft, 
dass dort vielmehr der Raketenantrieb mechanisch wesentlich iiberlegen ist. 

Wir bezeichnen mit m, die Fahrzeugmasse in einem beliebigen Augenblick, 
in dem das Fahrzeug die Geschwindigkeit v gegeniiber dem erdfesten Beobachter 
hat, und mit dm, die in diesem Augenblick dt verbrauchte kleine Treibstoffmasse 
des Heizwertes H, welch letzterer sich der Einfachheit halber hundertprozentig 
in Energie umsetzen mége, das heisst, der rein maschinelle, innere Wirkungsgrad 
wird zu Eins angenommen. 

_ Im Falle des radangetriebenen Fahrzeuges werde die Masse dm, ohne Relativ- 

geschwindigkeit zum Fahrzeug von Bord gegeben, im Falle des raketenange- 
triebenen Fahrzeuges jedoch mit der entgegengesetzt zur Fahrtrichtung wei- 
senden Relativgeschwindigkeit w = //2 g H. 

Die Fahrzeugmasse zu Anfang des Antriebvorganges sei m,, am Ende des An- 
triebes m,, der Treibstoffvorrat also m, — m,. 

Dann gilt fiir beide Falle: 


1. Radantrieb 
(Widerstands- und schwerefreie Beschleunigung) 


a) Erhaltung der Masse: 
dm, = —dm,; 


die abgegebene Masse dm, muss natiirlich gleich der Massenabnahme — dm, des 
Wagens sein. 
b) Erhaltung der Energie: 


My VU We gate, 
a( "2" ) = —dmy, ea + dm, = > 


die Anderung der kinetischen Energie d(m, v2/2) des Wagens muss gleich sein 
dem Energieverlust infolge der Abgabe der nach der Abtrennung noch mit der 
kinetischen Energie (dm, v?/2) ausgestatteten kleinen Masse dm, und dem Ener- 
giezuwachs (dm, w*/2) des ganzen Systems infolge der in diesem Augenblick frei 
gemachten Heizwertenergie (dm, g H) der verbrauchten Masse dmg. 

Durch Eliminieren von dm, aus beiden Gleichungen erhalt man die Differen- 
tialgleichung 

dm, 2 v dv 


elaine ae (1a) 


1) Arsenal de ]’Aéronautique, Chatillon-sous-Bagneux (France). 
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‘deren Integration die Grundgleichung des Radantriebes liefert :. 


Me Ze e-Vi2eH _ oe (1) 
Mg 
Der mechanische Antrieb P dt selbst folgt aus dem Impulssatz und Gleichung 
(la) zu: 
fs dv dm, w? gH dm, 
n ee ge 2G ee a (2) 


: 


das heisst, der Schub P an den Radern nimmt mit wachsender Fahrtgeschwin- 
‘digkeit hyperbolisch ab, die Antriebsleistung Pv ist konstant und unabhangig 
von der Fahrtgeschwindigkeit und bekanntlich nur durch Heizwert und sekund- 
lichen Treibstoffverbrauch bestimmt. 

Der dussere Antriebswirkungsgrad ist in jedem Augenblick: 


1 dm, v2 Mitt 0" I eae (3) 


A le eee 
dt 2 at 2 w2 


das heisst, der 4ussere Antriebswirkungsgrad des Radantriebes wird nur dann 
gleich Eins, wenn v?/2 g H gegen Null geht, was fiir unsere bekannten Radfahr- 
zeuge allerdings recht gut zutrifft. Dieser Antriebswirkungsgrad gilt auch fiir 
‘stationdre Fahrt, also gleichbleibende Fahrgeschwindigkeit, unter Uberwindung 
eines Fahrwiderstandes oder eines Schwerefeldes. 


2. Raketenantrieb 
(Widerstands- und schwerefreie Beschleunigung) 


a) Erhaltung der Masse: 
dm, = —dmy,. 
b) Erhaltung der Energie: ' 


2 2 2 
a( mas ) =—dm, 2 + dm, . 

Darin ist v, = w — v die Geschwindigkeit der abgestossenen Masse dm, nach 
der Abstossung relativ zum erdfesten Beobachter. 

Die Anderung der kinetischen Energie d(m,v?/2) des Raketenfahrzeuges muss 
gleich sein dem Energieverlust infolge Abgabe der nach der Abstossung noch 
mit der kinetischen Energie (dm, v3/2) ausgestatteten kleinen Masse dm, und 
dem Energiezuwachs (dm, w?/2) des ganzen Systems infolge der in diesem Augen- 
blick frei gemachten Heizwertenergie (dm, g H) der verbrauchten Masse dmy. 

' Purch Eliminieren von dm, und v, erhalt man die Differentialgleichung 


amy ee aw : (4a) 
my w 
deren Integration die bekannte klassische Raketengrundgleichung 
Me = o-V2eH = ea vlw (4) 
Mg 


liefert. 


te = 


, be by on Bn A ~*~ 
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Die Grundgleichung des Raketenantriebes wird gewéhnlich an Stelle des 
Energiesatzes aus dem anschaulicheren Impulssatz d(m, v) = dm, v, abgeleitet. 
Da der Impuls bekanntlich die Ableitung der kinetischen Energie nach der Ge- 
schwindigkeit ist, kommt man damit hier zum selben Resultat. ¥ 

Bemerkenswert ist bei diesen Ableitungen der Grundgleichungen, dass das 
totale Energiedifferential ; 

2 2 
a(m, =) = m,v dv + dm, 5 
im zweiten Glied den Energieverlust des Fahrzeuges infolge Abtrennung der mit 
der kinetischen Energie dm, v2/2 ausgestatteten Treibstoffmasse dm, = —dm, 
enthalt, und ferner, dass die absolute Grésse und der Verlauf der angewendeten 
Beschleunigungen ohne Einfluss auf das Ergebnis sind. 

Die Gegeniiberstellung der Gleichungen (1) und (4) in Figur 1 zeigt sofort, 
dass unter sonst gleichen Bedingungen zur Beschleunigung auf dieselbe Endge- 
schwindigkeit v beim Radantrieb wesentlich hGhere relative Treibstoffverbrauche 
(mz — ™,)/Mq, nétig sind als beim Raketenantrieb, sobald v/w> 1 wird. Bei ge- 
ringeren Endgeschwindigkeiten ist dagegen der Radantrieb giinstiger. Bei Werten 
von v/w, die gross gegen Eins sind, wird 


(m,/Mq) Rad ily —v2/w? 
(m,/™Mg) Rak 


das heisst, die Uberlegenheit des Raketenantriebes wachst rasch iiber alle Grenzen. 
Der mechanische Antrieb Pdi beim Raketenantrieb folgt aus dem schon er- 

wahnten Impulssatz und Gleichung (4a) zu 

dv dm, 


dm, .j=—> . 
Mm, FT ilins 7, V2eH, (5) 


iP 


das heisst, der Schub P der Rakete ist von der Fahrgeschwindigkeit unabhangig, 

die Antriebsleistung Pv nimmt proportional der Fahrgeschwindigkeit zu und 

kann grésser werden als g H dm,/dt, da als zusatzliche Energie noch die kine- 

tische Energie der an Bord mitgefiihrten Treibstoffmassen zur Verfiigung steht. 
Der aussere Antriebswirkungsgrad ist in jedem Augenblick bekanntlich: 


: iy 5 PY 2 v/w 6 
! dm, v? nt dm, w2 _ ‘whet } (6) 
dt 2 dive uae on 


Die zahlenmassige Gegeniiberstellung der 4usseren Wirkungsgrade nach den 
Gleichungen (3) und (6) des Radantriebes und des Raketenantriebes in Figur 2 
zeigt deren wesentliche Unterschiede, die auf die jeweils zweiten Glieder in den 
zitierten Energiesdtzen zuriickgehen: 

Der Energiezuwachs der Raketenfahrzeuge ist jeweils um 
dm 2 ) 

Sh (v8 = v8) = dm, 5 (= = 1) 

grosser als jener der Radfahrzeuge, das heisst, beide Vorgange sind gleichwertig 
bei v/w = 1/2; bei kleineren Fahrgeschwindigkeiten v ist der Raketenantrieb un- 
giinstiger als der Radantrieb, weil die kinetische Energie ihrer ausgestossenen 
Massen (relativ zum erdfesten Beobachter, Geschwindigkeitsrichtung entgegen 
der Fahrrichtung) grésser ist als jene des Radfahrzeuges. 
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